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SEPARATION  FROM  A  SLENDER  ROLLING  WING-BODY  COMBINATION 

by 

I.  P.  Jones* 

SUMMARY 

A  vortex-sheet  model  of  the  leading-edge  separation  is  used  to  study  the 
flow  about  a  slender  wing-body  combination  in  steady  rolling  motion  at  zero 
incidence.  The  equations  which  model  the  flow  are  derived  and  the  numerical 
method  to  solve  these  equations  is  described.  The  results  presented  concentrate 
on  those  values  of  the  parameters  for  which  comparison  may  be  made  with  the 
experimental  results  of  Harvey  for  a  rolling  delta  wing.  Results  obtained  using 
simplified  flow  models  are  also  described. 
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1  INTRODUCTION 

In  this  Report  we  consider  a  wing-body  combination,  consisting  of  a  slender 
wing  mounted  symmetrically  on  a  circular  cylinder,  and  rolling  steadily  about  the 
axis  of  the  cylinder,  which  is  aligned  with  a  uniform  stream.  The  flow  is  assumed 
to  separate  at  the  leading  edges  of  the  wing. 

When  a  plane  slender  wing  is  at  incidence  to  a  uniform  stream,  the  boundary 
layers  on  its  upper  and  lower  surfaces  leave  the  leading  edges  to  coalesce  into 
free  shear  layers  which  spiral  into  vortex  cores  lying  above  the  upper  surface  of 
the  wing.  The  flow  is  observed  to  be  steady  and  stable,  and  it  is  possible  to 
exploit  the  characteristics  of  the  leading-edge  vortices,  particularly  to  enhance 
the  lift  of  the  wing  at  low  speeds. 

If  the  wing  at  incidence  acquires  a  rolling  motion,  fast  enough  to  produce 
a  helix  angle  at  the  wing  tip  as  large  as  the  angle  of  incidence,  the  simplicity 
of  this  flow  pattern  is  destroyed.  At  some  station  along  the  ascending  leading 
edge  the  shear  layer  will  cease  to  feed  the  upper  surface  vortex  found  upstream 
and  will  instead  roll  up  into  a  second  vortex  below  the  wing.  Multiple  vortices 
will  also  form  from  the  leading  edges  of  a  wing  oscillating  in  roll  at  zero  inci¬ 
dence.  The  basic  simplicity  of  the  wing  at  incidence  is  preserved  in  two  cases: 
the  wing  at  incidence  and  with  small  disturbances  in  roll,  and  the  wing  at  zero 
incidence  in  steady  roll.  In  the  former  the  pair  of  Vortices  still  lie  above  the 
upper  surface  of  the  wing:  in  the  latter  they  lie  on  opposite  surfaces  of  the 
wing,  below  the  ascending  panel  and  above  the  descending  panel.  The  latter  prob¬ 
lem  was  chosen  for  the  present  investigation  because  the  nonlinear  effects  were 
expected  to  be  larger. 

The  conical  laterally  symmetric  problem,  associated  with  a  slender  delta 

wing  at  incidence,  has  attracted  considerable  attention  because  it  is  the  simplest 

problem  which  exhibits  the  characteristic  features  of  this  important  type  of  real 

1  2 

flow.  In  the  first  attempts  to  model  such  a  flow  ’  the  vortex  core  and  shear 
layer  were  replaced  by  an  isolated  potential  vortex  but  the  resulting  solutions, 
when  compared  with  experimental  results,  were  unsatisfactory.  To  overcome  the 
limitations  of  the  isolated  vortex  model.  Mangier  and  Smith^  included  some  of  the 
effects  of  the  shear  layer  by  incorporating  into  their  model  a  finite  vortex  sheet 
of  simple  form  which  extended  from  the  leading  edge  of  the  wing  towards  the  core 
of  the  vortex.  The  agreement  with  the  measured  lift  was  better  than  that  achieved 
with  the  earlier  models  but  the  vortex  position  was  still  not  predicted 
satisfactorily. 
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With  the  advent  of  high  speed  computers  Smith  obtained  the  first  satisfac¬ 
tory  solutions  to  the  original  problem.  The  model  he  used  was  similar  to  that  of 
Mangier  and  Smith,  except  that  the  shape  and  strength  of  the  vortex  sheet  were 
not  constrained  to  any  predefined  forms.  His  method  of  solution  reduced  the 
equations  which  express  the  boundary  conditions,  and  so  determine  the  position  and 
strengths  of  the  vortex  core  and  sheet,  to  a  set  of  nonlinear  algebraic  equations 
which  was  solved  using  a  nested  iteration  scheme.  The  numerical  scheme  owed  its 
origins  to  some  extent  to  physical  insight  and  associated  certain  variables  with 
certain  equations.  This  method  gave,  at  best,  linear  convergence.  For  low  values 
of  the  incidence  this  scheme  fails  to  converge  and  Barsby^,  using  a  more  accurate 
representation  of  the  sheet,  together  with  an  improved  iterative  procedure,  was 
able  to  extend  the  range  of  solutions  to  include  extremely  small  incidences. 

.  4 

The  assumptions  made  by  Smith  in  his  study  of  the  symmetric  conical  flow, 
which  are  described  in  more  detail  in  his  paper,  are  adopted  in  this  Report  and 
are  as  follows. 

(1)  The  flow  is  effectively  inviscid.  This  assumption  is  substantiated  by 
experimental  observations  which  indicate  that  viscosity  is  important  only  in  the 
thin  boundary  layers  on  the  surface  of  the  wing,  a  region  of  secondary  separation 
underneath  the  vortex  core,  in  the  shear  layers  and  in  a  small  region  of  the 
vortex  core.  The  experiments  confirm  that  most  of  the  large-scale  features  of 
the  flow  are  independent  of  the  Reynolds  number.  The  secondary  separation  which 
is  induced  by  an  adverse  pressure  gradient  caused  by  the  primary  vortex  is  ignored 
in  this  model. 

(2)  The  vorticity  in  the  fluid  is  condensed  onto  vortex  sheets  which  emanate 
from  the  leading  edges.  Each  sheet  is  infinite  in  length  and  rolls  up  into  a 
tight  spiral  core.  This  infinite  sheet  is  replaced  by  a  finite  outer  part 
springing  from  the  leading  edge  and  an  isolated  potential  vortex.  The  circulation 
on  the  spiral  core  is  condensed  onto  the  isolated  vortex  leaving  behind  a  trace. 
Across  this  trace  there  is  a  discontinuity  of  pressure  and  fluid  speed  that  is 
independent  of  its  shape.  This  trace  is  represented  by  a  cut  which  joins  the 
isolated  vortex  to  the  free  end  of  the  finite  vortex  sheet.  This  approximation 
appears  to  be  adequate  to  describe  the  flow  outside  the  core  region. 

20 

(3)  The  slender-body  theory  of  Munk,  Jones  and  Ward  is  used.  Thus  the 
velocity  potential  satisfies  Laplace's  equation  in  planes  perpendicular  to  the 
main  flow,  the  cross-flow  plane.  This  theory  neglects  the  influence  of  the 
trailing  edge  which  is  important  for  an  upstream  distance  of  the  order  of  the 


wing  span  at  low  speeds.  The  results  obtained  are  independent  of  Mach  number  but 
deteriorate  as  the  reduced  aspect  ratio  (A^l“  -  I)  increases  for  M  1  . 

The  effect  of  these  assumptions  is  to  reduce  the  problem  to  that  of  finding 
a  velocity  potential  which  satisfies  the  two-dimensional  form  of  Laplace's  equa¬ 
tion  subject  to  zero  normal  velocity  on  the  wing,  an  appropriate  condition  at 
infinity  due  to  the  incident  flow  and  certain  conditions  on  the  vortex  sheet,  the 
isolated  vortex  and  cut,  and  at  the  leading  edge. 

The  boundary  conditions  on  the  vortex  sheet  are  that  the  vortex  sheet  is  a 

stream  surface  of  the  three-dimensional  flow,  the  stream-surface  condition;  and 

that  the  pressure  is  continuous  across  the  vortex  sheet,  the  pressure  condition. 

It  is  impossible,  however,  to  satisfy  a  pressure  condition  on  the  isolated  vortex 

4  . 

and  cut  and  the  best  that  can  be  achieved  is  to  make  both  transverse  components 
of  the  overall  force  acting  on  this  system  zero.  Finally  a  Kutta  condition,  th.it 
the  velocity  is  finite  at  the  leading  edge,  is  applied. 

Using  this  vortex  sheet  model  several  authors*3  1 1  have  studied  a  number  of 
conical  flow  configurations. 

For  configurations  with  lateral  symmetry,  making  the  overall  force  on  each 

vortex  and  cut  vanish  automatically  ensures  that  there  is  no  net  rolling  moment 

9,10 

carried  by  the  fluid.  For  yawed  wings  it  is  possible  ’  to  choose  the  extent  of 
the  vortex  sheets  from  the  two  leading  edges  so  as  to  minimise  t  lie  net  rolling 
moment  in  the  fluid.  In  the  present  problem,  on  the  other  hand,  central  symmetry 
implies  that  the  rolling  moments  which  remain  on  each  vor lex-and-cut  system  are 
equal  and  cannot  combine  to  give  zero  overall  moment.  However  the  residual 
rolling  moment  in  the  fluid  is  small  compared  to  the  rolling  moment  on  the  wing 
and  is  ignored  in  the  present  treatment. 

More  fundamentally  the  problem  of  the  rolling  wing  is  not  conical,  even  for 
a  delta  wing  in  the  slender-body  approximation,  since  the  disturbance  produced  by 
elements  of  the  wing  depends  on  their  distance  trom  the  axis  ot  roll  and  this 
increases  as  the  distance  from  the  apex  increases  along  each  conical  ray.  This 

means  that  the  solution  in  each  cross-flow  plane  depends  in  a  more  complicated 

1  2  13 

manner  upon  the  upstream  conditions.  Karlv  attempts  ’  to  study  non-comcal 

'> 

flows  used  simplified  flow  models  based  on  that  ot  Brown  and  Michael  .  The 

boundary  conditions  on  a  non-conical  vortex-sheet  embedded  in  a  steady  potential 

14  ...  15 

flow  have  been  derived  by  Smith  .  Using  these  boundary  conditions  Smith  has 

obtained  similarity  solutions  for  two  classes  of  wing  for  which  conical  flow 

emerges  as  a  special  case.  More  recently  Clark'*  and  Jones"  have  developed 
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numerical  marching  schemes  for  dealing  with  slender  wings  of  arbitrary  planform. 

Other  models  have  been  developed  to  attack  non-conical  flow  problems,  for  example 
1  8 

by  Sacks  ct  al  ,  with  varying  degrees  of  success.  A  reasonably  accurate  method 

for  predicting  forces  and  moments  on  wings  with  edge  separation  is  based  on  the 

!  9 

leading-edge  suction  analogy  of  Polhamus  .  This  analogy  states  that  the  non¬ 
linear  component  of  the  force  acting  normal  to  the  wing  in  separated  flow  is  equal 
to  t he  suction  force  acting  in  Lhe  plane  of  the  wing  due  to  the  singular  behaviour 

of  the  attached  flow  at  the  sharp  edge.  Despite  the  success  of  the  analogy  it  has 

not  been  possible  so  far  to  provide  a  wholly  rational  basis  for  it. 

In  this  Report  we  are  primarily  concerned  with  the  separated  flow  over  a 

wing-body  combination  at  zero  incidence  rolling  steadily  about  its  axis  with 
angular  velocity  Si  in  a  free  stream  of  velocity  U  .  The  wing-body  combination 
consists  of  a  slender  delta  wing  of  semi-apex-angle  y  mounted  symmetrically  on  a 
circular  cylinder  of  radius  r  .  The  span  of  the  configuration  is  2s  .  This 
configuration  is  sketched  in  Fig  1.  The  mathematical  model  and  numerical  method 
are,  however,  applicable  to  a  slender  wing  of  general  planform  symmetrically 
mounted  on  a  circular  cylinder.  As  explained  later,  solutions  can  be  obtained 
which  approximate  those  for  isolated  wings. 

Relative  to  any  set  of  non-rotating  axes,  the  flow  considered  is  time- 
dependent.  However,  by  introducing  axes  which  rotate  with  the  coni' igurat  ion,  the 
problem  is  reduced  to  an  equivalent  problem  in  steady  flow  with  a  pro-determined 
distribution  of  rotation  in  the  fluid. 

In  section  2,  the  attached  flow  solution  to  this  problem  is  obtained,  using 
the  conformal  mapping  of  the  region  of  the  cross-flow  plane  occupied  by  the  tluid 
on  to  the  region  of  the  transformed  plane  exterior  to  a  circle.  The  separated 
flow  is  then  represented  by  introducing  vortex  sheets  and  isolated  vortices  into 
the  attached  flow  field,  together  with  their  images  in  the  circle  in  the  trans¬ 
formed  plane.  In  this  way,  the  boundary  conditions  on  the  wing  and  body,  satis¬ 
fied  by  the  attached  flow,  are  preserved.  Details  of  the  construction  of  the 
separated  flow  model  are  given  in  section  4. 

The  boundary  conditions  which  must  be  satisfied  on  the  non-conical  vortex 
sheet  in  the  rolling  motion  are  derived  in  section  3,  following  the  derivation  of 
Ref  14.  The  zero  force  condition,  expressing  the  vanishing  of  the  total  force  on 
the  isolated  vortex  and  cut  which  represents  the  core,  and  the  Rut  La  condition  at 
the  leading  edge  are  formulated  in  section  4.  These  conditions  take  the  form  of 
ordinary  differential  equations  which  arc  sufficient  to  define  the  rates  of 
change,  in  the  streamwise  direction,  of  the  position  and  strength  of  the  isolated 


vortex  and  of  the  shape  and  strength  of  the  vortex  sheets,  in  terms  of  the 
velocity  field  in  the  cross-flow  plane. 

Near  the  apex  of  the  exposed  wing,  point  A  in  Fig  1,  the  semi-span  of  the 
exposed  wing  is  small  in  comparison  with  the  body  radius,  so  that  the  body  side 
appears  as  an  infinite  reflection  plane.  Moreover,  the  relative  velocity  between 
the  wing  and  the  undisturbed  fluid  is  approximately  uniform  over  a  small  region, 
so  the  flow  near  A  resembles  that  of  a  delta  wing  at  an  angle  of  incidence  equal 
to  rsl/ll  .  Hence  the  flow  near  A  may  be  obtained  from  a  known  conical  solution. 
This  provides  the  initial  conditions  for  the  set  of  ordinary  differential  equa¬ 
tions  described  above,  so  that  the  solution  is  now,  in  principle,  determined. 

The  numerical  solution  of  these  equations  is  described  in  section  5.  The 
process  falls  into  a  number  of  distinct  stages.  The  representation  of  the  vortex 
sheet  in  the  cross-flow  plane  is  discretised,  so  that  the  unknowns  become  a  finite 
set  of  coordinates  and  strengths,  and  the  boundary  conditions  become  a  finite  set 
of  ordinary  differential  equations.  The  streamwise  derivatives  which  occur  in 
these  equations  are  represented  by  central  differences  in  an  implicit  formulation. 
This  formulation  requires  the  solution  of  a  set  of  simultaneous  nonlinear  equa¬ 
tions  for  the  values  of  the  coordinates  and  strengths  of  the  sheet  and  isolated 
vortex  in  each  cross-flow  plane,  each  set  of  equations  depending  on  the  solution 
already  found  for  cross-flow  planes  further  upstream.  These  simultaneous  non¬ 
linear  equations  are  solved  by  a  multi-dimensional  form  of  the  Newton-Raphson 
method.  Because  the  initial  conditions  obtained  from  the  conical  solution  for  the 
wing  at  incidence  are  not  exactly  appropriate  to  the  rolling  problem,  instabili¬ 
ties  can  arise  in  the  first  few  steps  of  the  downstream  integration.  These  are 
overcome  by  an  averaging  process. 

The  only  detailed  experimental  results  available  for  comparison  are  those  ol 
Harvey^*,  for  an  isolated  delta  wing  of  semi-apex  angle  10°.  He  obtained  measure¬ 
ments  in  three  cross-flow  planes  at  three  different  rates  of  roll,  in  an  essen¬ 
tially  incompressible  flow.  The  present  calculations  have  been  directed  towards 
obtaining  results  for  these  nine  cross-flows.  The  results  obtained  are  discussed 
in  section  6  in  relation  to  Harvey's  observations.  As  described  above,  the 
presence  of  the  body  is  an  essential  feature  of  the  calculation  method.  In 
principle,  the  body  can  be  made  as  small  as  we  please,  relative  to  the  overall 
wing  span  at  the  station  of  interest.  Tn  practice  this  means  that  the  conical 
starting  solution  corresponds  to  a  very  weak  separation,  for  which  difficulties 
arise”*  in  obtaining  the  appropriate  solution.  However  it  emerges  that  solutions 
virtually  independent  of  the  body  size  can  be  found  for  bodies  small  compared  with 
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the  local  span.  To  provide  an  indication  of  the  influence  of  body  size  on  the 

2 

separated  flow,  the  isolated  vortex  model  of  Brown  and  Michael “  was  also  formula¬ 
ted  for  the  wing-body  combination  and  results  calculated  for  comparison  with  the 

">2 

results  of  Hanin  and  Mishne  for  the  wing  alone.  Finally,  the  leading-edge 
suction  analogy  of  Polhamus  was  applied,  within  the  context  of  slender-body 
theory,  to  the  wing-body  combination.  The  result:  are  discussed  in  section  6. 

2  THE  COMPLEX  POTENTIAL  FOR  ATTACHED  FLOW 

Within  the  framework  of  slender-body  theory  we  are  concerned  with  solutions 
of  Laplace's  equation  in  the  cross-flow  plane.  Since  this  equation  is  linear  and 
the  boundary  conditions  on  the  wing  and  body  and  at  infinity  are  also  linear,  the 
velocity  field  of  the  separated  flow  may  be  constructed  by  superposing  a  represen¬ 
tation  of  the  leading-edge  vortices  on  the  attached  flow  solution.  In  this 
section  therefore  we  consider  the  attached  flow  problem  for  the  rolling  wing-body 
combination  described  in  the  last  section  and  illustrated  in  Fig  1. 

We  define  a  right-handed  set  of  Cartesian  axes  (x,  y,  z)  fixed  in  space  with 
the  origin  0  on  the  axis  of  the  body  in  the  cross-flow  plane  through  the  apex  of 
the  exposed  wing.  The  x-axis  is  along  the  axis  of  the  body  in  the  downstream 
direction  and,  at  the  instant  of  time  that  we  are  considering,  l  =  0  say,  the 
y-axis  is  in  the  direction  of  the  starboard  wing,  as  snown  in  Fig  1.  This  frame 
of  reference,  since  it  is  fixed  in  space,  we  shall  call  the  'inertial  frame'.  We 
also  find  it  convenient  to  work  with  another  frame  of  reference,  the  'rotational 
frame',  with  a  set  of  Cartesian  axes  fixed  in  the  body  such  that  at  the  same 
instant  of  time,  t  =  0  ,  these  two  sets  of  axes  are  coincident.  In  this  rota¬ 
tional  frame  the  flow  is  steady  but  it  is  not  a  potential  flow.  In  the  inertial 
frame,  however,  the  flow  is  an  unsteady  potential  flow  and  in  the  remainder  of 
this  section  we  shall  work  entirely  within  this  inertial  frame  at  the  instant 
t  =  0  . 

We  make  use  of  slender-body  theory  and  since,  in  the  absence  of  separation, 
no  vorticity  is  being  converted  downstream  the  velocity  field  in  each  cross-flow 
plane  is  entirely  determined  within  that  plane.  Hence,  for  the  flow  in  each 
cross-flow  plane,  we  may  define  a  complex  potential 

W,(Z)  =  *  +  iV 

A 

where  7,  =  y  +  iz  with  $  the  velocity  potential  and  the  stream  function. 

The  effect  of  the  free  stream  is  not  included  in  W,  . 

A 
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The  boundary  condition  which  must  be  satisfied  on  the  wing-body  combination 
is  that  the  normal  velocity  relative  to  the  wing-body  is  zero.  Hence,  on  the  wing 

—  =  Qy  for  z=0  and  r  ^  |y|  ^  s  (1) 

where  i?  is  the  angular  velocity  about  the  positive  x-axis.  Note  that  Fig  I  is 
drawn  for  a  negative  value  of  Si  .  On  the  body 

3$  „  *  fl  2 

—  =  0  for  /y  +  z  =  r  (2) 

3n 

where  n  is  the  unit  vector  normal  to  the  body.  At  infinity  the  fluid  is  at  rest 
and  so  we  have 

44  v  0  as  Z  -*•  ">  (3) 

dZ 

Using  various  methods,  several  authors  have  considered  problems  of  rolling 

23 

wing-body  and  cruciform-tail-body  combinations,  for  example,  Lomax  and  Heaslett  , 
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Adams  and  Dugan  and  Bobbit  and  Malvestuto  .  The  method  used  here  is  based  on 
that  presented  by  Adams  and  Dugan  for  a  rolling  cruciform-tail-body  combination 
and  introduces  a  distribution  of  doublets  on  the  surface  of  the  wing  to  satisfy 
the  boundary  conditions  (I)  and  (2). 

First  we  transform  the  region  of  the  cross-flow  plane  outside  the  cross 
section  of  the  wing-body  combination  onto  the  exterior  of  the  unit  circle  in  the 
transformed  plane.  By  a  Joukowski  transformation  the  cylindrical  body  can  be 
transformed  into  a  slit  lying  on  the  real  axis.  The  wing  itself,  since  it  lies  on 
the  real  axis  in  the  physical  plane,  is  merely  transformed  into  an  extension  of 
the  same  slit  along  the  real  axis  under  this  transformation.  Thus  the 
transformation 


Z*  =  -r-2 — 2  (Z  +  r'!/Z) 
s^  +  r~ 

is  such  that  the  wing-body  combination  is  represented  in  the  '/.*■  plane  by  the 
slit 


-  I  <  7*  I  . 


«r 


10 


!  T 


This  slit  is  mapped,  again  with  the  aid  of  a  Joukowski  transformation,  into  the 
unit  circle  | A |  =  !  .  The  complete  transformation  oi  the  wing-body  combination 
into  the  unit  circle  is  then  given  by 


where  A  is  the  point  in  the  final  transformed  plane  corresponding  to  the  point 
Z  in  the  physical  plane. 

In  this  transformed  plane  we  see  from  (A)  that  the  points  A  on  the  arc  of 
the  circle  representing  the  upper  surface  of  the  starboard  wing,  r  <  y  <  s  ,  are 
given  by 


where 


ie 


A  =  e 


8 ,  =  cos 


0  <  6  <  0 


2rs 


2.2 
s  +  r 


(5) 


Similarly  points  on  the  lower  surface  of  the  starboard  wing  are  given  by 


A  = 


-ie 

e 


and  on  the  upper  and  lower  surfaces  of  the  port  wing  by 


A 


-10 


and 


+10 


respectively,  for  0  <  9  <  6^  as  shown  in  Fig  2a. 

i  0 

Let  us  now  consider  a  source  of  strength  ra  at  a  point  A  =  e  , 

L  0 

0  <  8  <  tt  ,  on  the  unit  circle  and  a  sink  of  strength  m  at  the  point  A  =  e 
as  shown  in  Fig  2b.  The  complex  potential  for  this  source-sink  configuration  at 
some  point  Aq  in  the  transformed  plane  is 


lo  -  6 


ie 


co  e 


■10 


(6) 


From  this  we  can  easily  recover  the  classical  result  that  the  stream  function  is 
piece-wise  constant  on  the  unit  circle,  so  that  the  source-sink  combination 
induces  no  normal  velocity  on  the  contour  except  at  the  locations  of  the  source 
and  sink.  This  scurce-sink  pair  corresponds  to  a  doublet  of  strength  m  in  the 
physical  plane  where  we  again  have  zero  normal  velocity  everywhere  on  the  wing- 
body  except  at  the  location  of  the  doublet  itself.  Furthermore,  because  we 
effectively  have  a  source  on  the  upper  surface  of  the  wing  and  a  sink  on  the 
lower  surface  the  velocity  of  the  fluid  normal  to  the  wing-body  at  the  location  of 
the  doublet  is  in  the  same  direction  on  both  sides  of  the  wing,  as  is  required 
from  ( 1 ) . 

We  now  consider  a  distribution  of  doublets  along  the  wing.  This  automatic¬ 
ally  satisfies  the  condition  of  zero  normal  velocity  on  the  body  (2)  as  we  have 
seen,  and  the  boundary  condition  (3)  at  infinity.  From  the  continuity  of  mass 
flow  across  the  surface  we  require 

dm  =  2  dy  (7) 

3z 


where  dm  represents  the  strength  of  the  doublet  on  the  segment  dy  of  the  wing 
in  the  physical  plane  and  34>/3z  is  given  by  (1).  From  this  distribution  of 
doublets  on  the  wing  we  have,  at  the  point  ^  in  the  transformed  plane, 


W  (O  = 
A  0 


dm  log 


co  -  lw 


dm  log 


S  ~  C(y) 
f'0  "  ^(y) 


where  t,(y)  is  the  point  in  the  transformed  plane  corresponding  to  the  point  y 
on  the  wing  in  the  physical  plane  and  the  bar  denotes  the  complex  conjugate. 
Introducing  (1)  into  (7)  and  then  substituting  the  result  into  (8)  we  obtain  for 
the  attached  flow  potential 
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WA  = 

I - - 1 

•K 

O 

00 

2  2  1 
co  "  c 

•c 

II 

Cfi 

61 
+  f 

fixil 

i  "  2Ue10  -  2Ue"10  | 

2  -2 
?o  -  5  J 

u 

II 

>> 

J 

0 

2n  j 

2 

C0  ~ 

2  2-2 

C  "  5  1 

since,  from  (5) ,  0  = 

0  when 

y  = 

s  and  6 

=  91 

when  y  =  r  . 

From  (4)  and  (5) 

,  since 

z  = 

y  on 

the 

wing. 

we  have 

2 

2yr 

cos  0 

2 

0  . 

y  _ 

cos  0  , 

“  +  Y 

The  integrals  in  (10)  may  now,  with  the  aid  of  (11),  be  evaluated,  using  standard 
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techniques.  The  details  of  the  integration  may  be  found  in  Jones  .  The  complete 
expression  for  W  is  given  by 

A 


2  /  2  2l6l 

nr  /  £  -  e 

TT  lo*  - - r2le' 


ifir2(c^ 


D9, 


2 

.4  ~  e 


2  2 
2tt£  cos  0, 


or2  9  /(4 

- 5 —  (Jx  +  sin  0  ]  log 


IT  COS  0 


1 )  +  i  tan  0  j ( C  +  1 ) 
(C2  -  O  -  i  tan  6] (C2  +  I), 


> 


4  2 

inu  -  I)r 

2  2 
4;  cos  0. 


2  2 
lflr  Q 


4  cos  0 


2  cos  20  ^  +1 


(12) 


where  x  =  4  +  1/4^  and  the  suffix  zero  has  been  dropped  from  . 

By  differentiating  this  expression  for  constant  x  ,  the  corresponding  com¬ 
plex  velocity  in  the  transformed  plane  may  easily  be  found  as 


dW 
_ A 

dr. 


2h?r'  sin  20 


tiir  0 


'**"  I  . I 

7~i  i  ’  \  +  2 

nf,  I  +  — —  -  2  cos  20  J  r.  ii  cos  0 


A 


^  7  ('•’  '  J5)  - 

COS  0 .  \  i,  / 


U'  -  I )  +  i  tan  0  (  (r.'  +  1 ) 
U'  -  1)  -  i  tan  0 ,  (?,'  +  1) 


4hlr'  tan  0 


u  cos  0 , 


l(l/2  1  \  .  2,  |  f. 

-  *•  +  —  I  +  *l"  }  7— > - ^ - 3 - 1 - 7 

('  \  r.-/  )  (U‘  -  D'  +  tan'O  (r'  +  O' 


h:r'  /  2  I  \  i i:r '  / 

- — -  U  +—  + - ~ ~ 3  A 

2  cos  0,i'.  '  r,  /  2  cos  O  r, 


-  2  cos  20  1 r  +1 


)  •) 

is’r  '  (,\'  +  1) 


(i,“  -  cos  20.) 


2  cos  0 . 


/  r."'  -  2  cos  20^,'  +  1 


J 

■  03) 


From  this  comp l i cat  oil  expression  (13)  it  may  be  verified  directly  that  we 
have  the  required  complex  vel  00  i  t  y  i  11  l  hi*  transformed  plane  which  satisfies  condi¬ 
tions  (I)  on  the  wing,  (2)  on  the  body  and  (3)  at  infinity.  Care  must,  however, 
he  taken  when  evaluating  t  lie  complex  logarithms  and  square  roots  in  (13)  to 
accommodate  the  changes  front  one  branch  to  another  which  may  be  encountered  011 
that  part  of  the  boundary  in  the  r.-plane  which  corresponds  to  the  wing.  This  is 
because  for  a  point  r.  on  the  boundary  the  integrand  in  (9)  becomes  singular. 
Hence  to  calculate  3W^/3r,  at  r.  =  eUl  ,  say,  from  (13)  we  must  evaluate 


1  im 
1-  -0 


dW 
_ A 

d  f, 


(O 


>e1(') 


where  t  is  real  and  positive,  to  ensure  that  the  correct  values  are  obtained. 

Knowing  the  attached  flow  potential  and  the  corresponding  complex  velocity 
through  (12)  and  (13)  we  are  now,  as  mentioned  in  section  1,  in  a  position  to 
construct  the  complex  velocity  in  the  presence  of  loading-edge  separation.  This 
is  carried  out  in  section  4.  First,  however,  in  section  3  we  derive  the  boundary 
conditions  to  be  satisfied  on  a  steadily  rolling  vortex  sheet. 


K 
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3  THE  BOUNDARY  CONDITIONS  ON  A  STEADILY  ROLLING  VORTEX  SHEET 

In  this  section  we  discuss  in  detail  the  boundary  conditions  to  be  satisfied 
on  a  vortex  sheet  which  is  rolling  steadily  about  an  axis  parallel  to  an  oncoming 
stream.  As  we  have  discussed  in  section  1,  the  conditions  that  must  be  satisfied 
are  continuity  of  pressure  across  the  sheet,  t  lie  pressure  condition,  and  the  con¬ 
dition  that  the  vortex  sheet  is  a  stream  surface  of  the  three-dimensional  flow, 
the  stream-surface  condition. 

In  a  rotating  frame  of  reference  the  motion  is  not  a  potential  flow.  How¬ 
ever  if  we  consider  a  set  of  axes  (x,  y,  z)  in  the  inertial  frame  defined  in  the 
previous  section,  that  is  fixed  in  space,  with  the  wing  rotating  about  the  x-axis 
with  constant  angular  velocity  if  =  (11,  0,  0)  ,  then  relative  to  these  axes  the 
flow  is  an  unsteady  potential  flow.  Hence,  from  Bernoulli’s  theorem  we  have 


3 1 


+ 


Ju 


U  +  £  =  c  (  L ) 

P 


U  4) 


where  u  =  V<t>  .  Thus,  from  continuity  of  pressure  across  the  sheet  we  have 


A(lf  +  *-  '  = 


(  I  ri) 


where  A  is  the  difference  operator  across  the  seoet,  'inside  minus  outside'. 

In  order  to  relate  3<t>/3t  to  the  velocity  field,  we  first  introduce 
Cartesian  axes  in  the  rotating  frame,  fixed  to  the  wing,  and  therefore  rotating 
about  the  x-axis  of  the  inertial  frame  with  constant  angular  velocity  s l  .  If  the 
new  y'  and  z'  axes  coincide  with  the  inertial  y  and  x  axes  at  time  t  =  0 
the  coordinates  are  related  by: 


y'  =  y  cos  Sit  +  z  sin  Sit 
z'  =  z  cos  Sit  -  y  sin  Sit 
y  =  y'  cos  Sit  -  z'  sin  sit 
z  =  y*  sin  sit  +  z'  cos  Sit 


(lb) 


The  velocity  vector  q  ,  relative  to  the  new  axes  is  given  by 

u  =  VI1  =  q  *  ll  a  r  (17) 


L 


i 


039 
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(Greenspan)  ,  where  r  is  the  radius  vector  from  the  origin  and  the  gradient 
operator  is  relative  to  the  inertial  frame.  Since  the  flow  is  steady  in  the  new 
axes,  we  write 


so  that 


<I>(x,  y,  z,  t)  =  F(x ,  y',  z') 


3$  =  3F_  3^1  3£_ 

at  "  ay'  at  a?.'  at 


=  si 


/  ,  3F  ,  DF  \ 

\z  w  " y 


by  (16). 

Similarly  we  find 


(18) 


and 


Hence,  by  (18), 


3F  _  3y_  a_i  3z 

3y'  3y  3y'  +  3z  3y' 


,  34>  ,  3-1 

=  cos  sit  —  +  sm  Sit  — 
3y  3z 


3F  .  3*  ^  34> 

3^  =  '  81n  3^  +  COS  SH  3i 


ai> 

at 


/  a*  a*  \ 

S:  (Z  3y  Y  3z  ) 


(19) 


Since  S2  =  (Si,  0,  0)  ,  equation  (17)  can  be  written  as 


q  =  u  -  (0,  -Slz,  sly),  with  u  =  V>1 


Now,  forming  the  scalar  product  and  taking  the  difference  across  the  sheet  at  the 
same  point  (y,  z) ,  we  have 


A(q  .  q)  =  A 


(s  •** 2i:z  #  '  £) 
A(a  •  a  *  2  If) 


=  0 


by  ( I  9) 
by  (IS) 


We  may  therefore  replace  the  pressure  condition  (15)  by 
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A(q  .  q)  =0  .  (20) 

In  other  words  the  jump  in  the  magnitude  of  the  velocity  in  the  rotational  frame, 
in  which  we  find  it  convenient  to  work,  is  zero. 

The  second  boundary  condition,  as  mentioned  earlier,  is  that  the  vortex 
sheet  is  a  stream  surface.  Thus 

q  .  n  =  0  (21) 

where  n  is  the  unit  vector  normal  to  the  surface. 
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These  conditions  (20)  and  (21)  are  now  similar  to  those  used  by  Smith  in 
his  derivation  of  the  boundary  condition  for  steady  potential  flow,  except  that 
the  actual  velocity  u  which  appears  in  his  conditions  is  replaced  by  the  rela¬ 
tive  velocity  q  in  the  rotating  frame.  For  this  reason  the  following  analysis 
follows  closely  that  of  Smith  although  some  of  the  detail  is  omitted.  Our  aim  is 
to  express  the  boundary  conditions  (20)  and  (21)  in  forms  which  are  suitable  for 
a  numerical  calculation.  The  first  step  is  to  express  them  in  terms  of  quan¬ 
tities  defined  on  the  section  of  the  vortex  sheet  by  a  plane  normal  to  the  free 
stream. 

We  define  a  coordinate  system  (£,  q,  40  with  £  as  a  coordinate  measured 
in  the  direction  of  the  uniform  stream,  and  with  q  and  \Jj  as  coordinates  in 
th°  cross-flow  plane  normal  to  the  uniform  stream.  The  vortex  sheet  T.  ,  being  a 
stream  surface,  cannot  lie  entirely  in  a  plane  £  =  constant  and  the  coordinate 
system  is  chosen  so  that  I  does  not  lie  in  a  surface  q  =  constant  .  This  con¬ 
figuration  is  shown  in  Fig  3.  We  define  a  and  b  as  unit  vectors  in  the 
tangent  plane  to  the  sheet  with  a  in  the  direction  of  constant  q  and  increas¬ 
ing  £  ,  and  b  in  the  direction  of  constant  £  and  increasing  q  .  At  the 
point  P  which  we  are  considering  a  ~  b  is  normal  to  the  sheet  and  hence  the 
stream  surface  condition  (21)  q  .  n  =  0  ,  may  be  written  as 

[a  a  b|  .  q |  =  I  a  a  b  ]  .  q^  =  0  (22) 

where  q  is  the  velocity  vector  relative  to  the  rotating  frame  and  the  subscript 
I  refers  to  quantities  on  the  inside  of  the  sheet  with  reference  to  the  direction 
of  the  vector  a  /s  b  ,  and  the  subscript  2  to  quantities  on  the  outside  of  the 
sheet . 

We  can  form  a  set  of  orthogonal  axes  with  unit  vectors  £  ,  b  and  i  ^  h  , 
where  i  is  the  unit  vector  in  the  strcamwise  direction  and  consequently  j  „  b 


fstn 


U39 


T 


I  7 

is  directed  towards  the  inside  of  Y.  along  the  normal  to  the  section  of  Y.  by 
the  plane  f  =  constant  .  If,  therefore,  we  write  q  as 

q  =  (U  +  q i  +  qtb  +  q  i  /v  b  (23) 

and  take  the  scalar  product  of  this  equation  with  a  b  ,  we  obtain  from  (22) 

(a  ~  b)  .  i(U  +  q^.)  +  (a  .  i)qn  =  0  (24) 

since  (a  *  b)  .  b  =  0  and  (a  b)  .  (i  a  b)  =  a  .  i  .  For  small  disturbances 
q  «  U  we  obtain  for  the  stream  surface  condition  (24): 


(a  b)  .  i 


U 


a  .  i 


(25) 


a  condition  which  applies  on  both  sides  of  the  vortex  sheet. 

The  pressure  condition  (20)  may  now,  from  (23),  be  written  as 


A 


<  a  2  2  1 

{(U  +  V  +  qt  +  qn}  = 


=  0 


We  have  then,  again  assuming  q  «  U  , 


UAqf,  +  qtmAqt  =  ° 


(26) 


where  q  is  the  mean  velocity  component  tangential  to  the  section  of  the  sheet 
by  the  cross-flow  plane: 


tm 


=  Hqu 


qt2> 


and  the  continuity  of  q^  across  the  sheet,  expressed  by  (25)  has  been  used. 

Since  the  component  of  relative  velocity  normal  to  the  sheet  is  zero,  the 
velocity  difference  across  the  sheet,  Aq  ,  lies  in  the  tangent  plane  to  the 
sheet  and  hence  can  be  written  as 


Aq 


\a  +  )ib 


!  8 


where  X  and  u  are  scalar  quantities.  Taking  the  scalar  product  of  this  equa¬ 
tion  with  a  and  b  in  turn  we  may  calculate  X  and  we  find  that 

x(l  -  (a  .  b)“)  =  a  .  Aq  -  (a  .  b)(b  .  Aq) 

Hence,  since  i  .  b  =  0 

Aqr  =  i  .  Aq  =  Xi  .  a 

=  (a  .  Aq  -  (a  .  b)(b  .  Aq) ) ( i  .  a)/(l  -  (a  .  b)2)  .  (27) 


Now,  from  (17) 


a 


2 1 


a  .  (U|  -  ij  „  r) 


h.  dr. 


a  .  (iJ  ^  r) 


(28) 


where  h  d(  is  the  arc  length  along  the  curve  rj  =  constant  on  X  and  the 

partial  derivative  is  for  constant  n  on  X  .  This  expression  differs  from  that 
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found  by  Smith  for  the  steady  potential  flow  case  by  the  term  which  involves 
the  vector  S2  explicitly.  However,  when  we  take  differences  across  the  vortex 
sheet  these  terms,  being  functions  of  position  only,  disappear  and  we  have,  from 


(28) 


a  .  Aq 


\  _  J_  JIM 

3f,  j  hf.  3f, 


(24) 


since  the  differentiation  3/3f  is  along  a  curve  lying  in  the  vortex  sheet  X  . 
Here  A<t>  is  the  circulation  or  the  jump  in  potential  across  the  vortex  sheet. 

In  a  similar  manner,  from  (23) 


Aqt 


b  .  Aq 


_l_  3M 

h  3r| 
n 


(30) 


where  h^dn  is  the  arc  length  along  the  curve  f;  =  constant  on  X  .  Substitu¬ 
ting  (29)  and  (30)  into  (27)  and  using  the  resulting  equation  for  Aqf  in  (26), 
we  obtain  for  the  pressure  condition: 


_l_  3M> 

hf  3f. 
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(31) 
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The  two  conditions  (25)  and  (31)  are  the  same  as  those  derived  by  Smith 
except  that  the  actual  velocities  are  replaced  by  the  velocities  relative  to  the 
rotating  frame.  In  his  paper  Smith  went  on  to  express  these  conditions  for  a  set 
of  polar  coordinates  but  in  the  present  work  we  choose  to  express  .nem  in  a 
different  set  of  coordinates.  Barsby^  has  shown  that  an  intrinsic  coordinate 
system  is  better  able  to  handle  the  shapes  of  leading-edge  vortex  sheets  that  are 
encountered  at  low  incidences.  For  this  reason  the  coordinates  n  and 
are  chosen  in  the  following  way.  f,  is  identified  with  x  ,  the  coordinate 
measured  in  the  downstream  direction,  so  that  h  =  1  .  >)  is  defined  by 


where  o  is  the  arc  length  along  the  sheet  measured  in  the  plane  x  =  constant  , 
and  C  is  an  arbitrary  function  of  x  .  We  should  like  to  keep  the  angle  in  the 
cross-flow  plane  which  the  line  joining  the  isolated  vortex  to  the  free  end  of 
the  vortex  sheet  makes  with  the  wing  approximately  constant,  and  to  generate  this 
vortex  sheet  by  a  fixed  range  of  values  of  t)  .  To  achieve  this  we  choose  C  to 
increase  with  x  to  keep  pace  with  the  growing  size  of  the  vortex  system,  if  is 
the  angle  made  by  the  tangent  to  the  vortex  sheet  with  some,  as  yet  unspecified, 

reference  line  in  the  cross-flow  plane.  We  also  define  a  set  of  Cartesian  axes 

(x,  y,  z)  with  y  in  the  direction  of  if'  =  0  and  with  unit  vectors  |  j  and 

k  in  the  direction  of  the  axes. 

In  this  case,  in  the  expressions  for  the  stream  surface  and  pressure  condi¬ 
tions  (25)  and  (31)  i  is  as  we  have  just  defined  it,  a  is  a  unit  vector  in 
the  tangent  plane  to  in  the  direction  of  n  constant  and  x  increasing,  and 

b  is  a  unit  vector  in  the  tangent  plane  with  i p  constant  and  n  increasing. 
These  are  sketched  in  Fig  3c.  Thus 

(33) 

1  14) 


where  3/3x1  denotes  the  downstream  derivative  on  the  surface  !'  ,  keeping  n 
n 

constant.  By  the  s 1 ender-body  assumptions: 
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and  hence,  from  (34): 
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Since  o  is  the  arc  length  in  the  direction  of  n  increasing,  h  dti  =  do  and 


so,  by  (30); 
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Using  this,  and  equations  (33)  and  (35),  we  obtain  from  (25)  and  (31)  the  stream 
surface  and  pressure  conditions  in  the  forms: 
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These  equations  will  also  be  applicable  to  the  steady  potential  flow  case, 

in  which  the  velocities  relative  to  the  sheet  in  the  cross-flow  plane,  q  and 

tm 

q  ,  become  the  actual  velocities  (3T/3o)  and  3T/3n  .  Thus  in  this  intrinsic 
n  m 

coordinate  system  the  relevant  boundary  conditions  for  a  steady  potential  flow, 
lor  example  a  wing  at  incidence,  are 
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For  conical  flow,  if  we  choose  >  t<>  he  the  wing  semi-span  so  that  the  vortex 

sheet  may  grow  linearly,  we  have 
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and 


tan  "y 
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8A$ 

9x 


z 

tan  y  — 
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M 

tan  y  — 
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where  s  is  the  semi-span  of  the  delta  wing  and  y  is  the  wing  semi-apex  angle. 
Hence,  the  conditions  (38)  and  (39)  become 


3*  U  tan  y  <  ,  \  /  /  r\\ 

—  - - - - -  (y  sin  4>  -  z  cos  if)  (40) 

and 

A*  “  lr  ( y  cos  *  +  z  sin  *  ■  u~rlr7  (B)n)  •  (41) 

If  we  finally  consider  a  polar  coordinate  representation  of  the  vortex  sheet  in 
the  cross-flow  plane,  we  have 


and 


r  cos  <(>  =  y  cos  if  +  z  sin  if 


r  sin  if  =  y  sin  if  -  z  cos  4* 


where  if  is  the  angle  made  by  the  radius  vector  with  the  tangent  to  the  sheet. 
Then,  from  (40)  and  (41)  we  have 


and 


U  tan  y  r  sin 


.  .  9A4>  . 

M  =  — —  I  r  cos  * 

9o  * 


(r  cos  * '  (i)J  ■ 


which  are  the  usual  conditions  for  a  conical  flow  . 

4  THE  EQUATIONS  GOVERNING  THE  SEPARATED  FLOW 

In  the  two  preceding  sections  we  have  obtained  an  expression  for  the  com¬ 
plex  potential  of  the  attached  flow  in  section  2  and  derived  convenient  forms  for 


i. 
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the  boundary  conditions  to  be  applied  on  the  vortex  sheet  in  section  3.  In  this 
section  we  add  terms  representing  the  leading-edge  vortices  to  the  attached  flow 
potential  to  obtain  a  complete  description  of  the  separated  flow  field;  and  we 
formulate  the  remaining  conditions  which  this  flow  field  must  satisfy,  the  Kutta 
condition  at  the  leading  edge  and  the  vanishing  of  the  force  acting  on  the 
coordinates  of  the  isolated  vortex  and  cut.  Finally  the  complete  system  of 
equations  is  presented  in  a  non-dimensional  form  based  upon  the  radius  r  of  the 
cylindrical  body  and  the  free  stream  speed  U  . 

4 .  1  Complex  potential  for  separated  flow 

As  described  in  the  Introduction  the  fully  rolled-up  vortex  sheet  of 
infinite  extent  which  is  sketched  below  the  port  wing  in  Fig  1  is  replaced  by  a 
finite  length  of  vortex  sheet,  an  isolated  vortex  filament  and  a  cut  connecting 
these,  as  illustrated  above  the  starboard  wing  in  Fig  I.  The  configuration  to  be 
represented  has  central  symmetry  about  the  axis  of  the  cylinder,  so  that  it  is 
sufficient  to  specify  the  shape  and  strength  of  the  vortex  system  in  the  first 
quadrant . 

In  section  2  we  found  expressions  for  the  complex  potential  (12)  and  its 
derivative  (13)  at  time  t  =  0  which  satisfy  the  boundary  conditions  on  the 
cross-section  of  wing-body  combination  and  at  infinity  and  which  are  analytic  on 
the  cross-flow  plane  outside  the  cross-section.  We  now  wisli  to  add  terms  which 
have  an  appropriate  singular  behaviour  on  the  sections  of  the  vortex  sheets  and 
at  the  positions  of  the  vortex  filaments  and  which  do  not  disturb  the  boundary 
conditions  on  the  wing-body  combination  and  at  infinity.  This  is  easily  achieved 
by  using  the  conformal  transformation  (4)  of  the  cross-flow  plane  of  Fig  1  onto 
the  4-plane  of  Fig  2.  The  wing-body  contour  becomes  the  unit  circle  in  the 
4-plane  and  so  the  boundary  condition  on  the  contour  is  preserved  by  including 
the  imago  inside  the  circle  of  every  singularity  introduced  outside  it.  The 
boundary  condition  is  also  preserved  if  a  vortex  of  arbitrary  strength  is  added 
at  the  centre  of  the  circle  in  the  4-plane.  However,  since  the  total  circulation 
of  each  element  of  vorticity  in  the  flow  field  and  its  image  inside  the  circle  is 
zero,  a  vortex  at  the  centre  of  the  circle  represents  a  net  circulation  about  the 
cross-sec t ion  of  the  configuration.  No  such  circulation  can  exist,  since  any 
contour  surrounding  the  cross-section  can  be  displaced  upstream,  following 
streamlines,  into  the  undisturbed  flow.  Hence  no  vortex  should  exist  at  the 
centre  of  the  circle. 

The  required  velocity  in  the  cross-flow  plane  is  then  the  complex  conjugate 
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where  aw^/34  is  given  by  (13).  r  is  the  circulation  of  the  isolated  vortex 
filament  in  the  first  quadrant,  this  vortex  lying  at  the  point  4^  in  the 
4-plane.  The  arc  length  along  the  vortex  sheet  in  the  first  quadrant  in  the 
physical  plane,  o  ,  is  used  as  the  integration  variable  and  takes  the  value  o 

E 

at  the  free  end  of  the  sheet,  A4>  is  the  jump  in  potential  across  the  sheet  and 
the  minus  sign  arises  because  A$  decreases  as  a  increases,  and  4(0)  is  the 
position  of  a  current  point  of  the  starboard  sheet  on  the  4-plane.  The  success¬ 
ive  terms  in  each  curly  bracket  in  (42)  arise  from  the  starboard  singularity,  its 
image  in  the  unit  circle,  the  port  singularity  and  its  image  in  the  circle. 

4.2  Conditions  on  the  separated  flow  field 

The  construction  of  relation  (42)  ensures  that  the  boundary  conditions  on 
the  wing-body  and  at  infinity  are  satisfied.  The  boundary  condition  on  the  vor¬ 
tex  sheets  (and  as  pointed  out  above,  it  is  enough  to  consider  one  sheet)  are 
given  by  equations  (36)  and  (37),  in  terms  of  the  velocity  q  relative  to  the 
rotating  frame  of  reference  in  which  the  flow  is  steady,  but  rotational,  q  is 
related  to  the  irro^ational  velocity  field  u  of  the  unsteady  potential  flow  at 
time  t  =  0  by  equation  (17).  Taking  components  of  (17)  along  the  tangent  and 
normal  to  the  cross-section  of  the  sheet  we  have 
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Combining  these  two  into  a  single  complex  equation  we  have 
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Integrating  this  last  equation  gives  a  convenient  expression  for  the 
coordinates  of  the  sheet  in  the  physical  plane: 


n 

Z  =  s(x)  +  J  e1 ’ do 
0 


(44) 


The  Kutta  condition  at  the  leading  edge  requires  that  the  velocity  there  is 
finite.  Since  in  the  transformation  (4)  the  leading  edge,  Z  =  s  ,  at  which 
r,  =  1  ,  is  a  singular  point  at  which  3 Z / 3 c  vanishes,  the  Kutta  condition 
implies  that 

■—  =  0  at  r,  =  1  .  (4  5) 

O  L, 


We  now  consider  the  forces  acting  on  the  cut  and  on  the  isolated  vortex 
filament.  The  difference  in  pressure  across  the  cut,  inside  minus  outside,  is 
given  by  (14)  as 


Ap 


Equation  (19)  expresses  3$/3t  in  terms  of  the  cross-flow  velocities  and  these 

■> 

are  continuous  across  the  cut.  Hence,  neglecting  (3v/9x)~,  the  pressure  jui 
becomes 


Ap 


=  -  pUA 


9x 


since  the  jump  in  potential  across  the  cut  is  just  the  circulation  about  the  vor¬ 
tex  filament.  The  force  on  the  section  of  the  cut  between  two  cross-flow  planes 
a  small  distance  Ax  apart  produced  by  this  pressure  difference  is  given  in 
magnitude  and  direction  by 


-  i(Zy  -  ZE)  ApAx  =  iPU(Zv  -  ZE)  AX 


(4b) 


where  and  Z^  are  the  position  of  the  vortex  filament  and  the  free  end  of 

the  vortex  sheet. 

Randall'^  has  shown  that  the  force  on  an  element  of  a  vortex  filament  of 
length  Ax  in  an  unsteady  potential  flow  is  given  by  the  classical  Joukowski 
formula  for  steady  flow,  viz: 
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-  ipFQAx 

where  Q  is  the  complex  velocity  normal  to  the  vortex  element.  In  the  present 
situation  there  are  three  contributions  to  Q  .  One  comes  from  the  unsteady 
potential  fLow  in  the  cross-flow  plane  and  is  obtained  by  deducting  the  influence 
of  the  vortex  itself  from  the  complete  cross-flow  velocity,  giving 


2ni  (Z  -  Zv) 


where  the  bar  denotes  the  complex  conjugate.  The  second  contribution  comes  from 
the  relative  velocity  between  the  rotating  and  inertial  frames,  given  from  (17)  as 


iOZy  . 


The  third  contribution  comes  the  inclination  of  the  vortex  filament  to  the  free 


stream,  viz: 


9,  =  -  u  -nr 


dx 


Writing  Q  =  Qj  +  ^2  +  -3  ’  anci  e9uat^n8  to  zero  the  sum  of  the  force  on  the 
element  of  the  vortex  and  force  (46)  on  the  element  of  the  cut,  we  obtain 
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This  is  the  final  condition  to  be  imposed  on  the  vortex  system. 

4 . 3  Non-dimensional  form  of  the  equations 

We  define  a  non-dimensional  parameter  P  to  describe  the  rate  of  roll 


P  =  - 


The  minus  sign  is  introduced  so  that  we  have  a  vortex  sheet  in  the  first  quadrant, 
as  in  Fig  1,  where  P  is  positive.  If  the  configuration  is  regarded  as  moving 
forward  with  speed  U  ,  the  apex  A  of  the  exposed  wing  describes  a  helix  whose 
inclination  to  its  axis  is  tan  'p  . 
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Since  the  body  radius  r  is  independent  of  x  ,  it  is  convenient 
as  a  length  scale,  with  the  free  stream  velocity  U  as  a  velocity  scale 
this  point  on,  the  same  symbols  as  before  will  be  used ,  but  they  will  be 
the  quantities  made  dimensionless  in  this  way. 

In  non-dimensional  form  then  the  stream-surface  condition  (36)  and 
condition  (37)  become 
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where  now,  from  (43)  and  (48) 
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The  force  condition  (47)  becomes 
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and  the  Kutta  condition  (45)  remains  as 
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The  complex  conjugate  of  the  attached-flow  velocity  (13)  is  now  given  by 
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and  the  complex  conjugate  of  the  velocity  of  the  separated  flow  (42)  is 
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(54) 


The  formula  for  the  sheet  coordinates  (44)  remains  unaltered  in  non-dimensional 
form  and  so  the  relationship  between  (y,  z)  and  (iji,  o)  may  be  written,  as  before. 


Z  =  s(x)  +  j  e‘’do 
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With  the  complex  velocity  in  the  transformed  cross-flow  plane  now  given  by 
(,54)  the  equations  which  express  the  boundary  conditions,  namely  (49),  (50),  (52) 
and  (55),  and  so  determine  the  position  and  strengths  of  the  isolated  vortices 
and  vortex  sheets,  may  be  reduced  to  a  set  of  nonlinear  ordinary  differential 
equations  in  the  streamwise  variable  x  .  Since  the  flow  near  the  apex  of  each 
exposed  wing  panel  is  approximately  conical,  initial  values  of  the  unknowns  are 
available  from  previous  solutions,  and  the  equations  can  be  solved  using  a 
marching  procedure.  This  procedure  is  described  in  the  following  section. 

For  a  given  combination  of  a  slender  wing  on  a  cylindrical  body,  ;V  for  a 
specific  non-dimensional  local  semi-span  s(x)  ,  the  solution  depends  only  on  the 
parameter  P  defining  the  rate  of  roll. 

5  FINITE  DIFFERENCE  FORMULATION  AND  SOLUTION  PROCEDURE 

The  finite  difference  formulation  that  we  use  to  solve  the  nonlinear 
differential  equations  derived  from  expressing  the  boundary  conditions  as  in 
section  4  can  be  split  into  two  separate  parts.  The  first  is  the  discretization 
in  the  streamwise  direction.  The  second  concerns  the  numerical  representation  of 
the  vortex  system  in  each  cross-flow  plane  and  Lhe  solution  of  the  resulting 
equations  for  each  cross-flow  plane. 

In  the  finite  difference  representation  that  we  adopt  the  shape  of  the 
sheet  in  each  cross-flow  plane  is  defined  by  the  value  of  the  angle  C  at  n 
points  on  it,  with  the  siieet  strength  defined  by  -  3A<f/;)o  ,  again  evaluated  at  the 
same  points.  These  points  are  mid-points  of  specially  chosen  intervals  and  for 
this  reason  they  are  known  as  intermediate  points.  The  intervals  are  defined 
later.  The  values  which  represent  the  strength  and  shape  of  each  sheet,  together 
with  the  real  and  imaginary  parts  of  the  position  Z^  ,  and  t  tie  strength  ;  ,  of 

the  corresponding  isolated  vortex  give  2n  +  3  unknowns  represented  by  the 
vector  X  ,  say,  to  be  determined  in  each  cross-flow  plane. 

Near  its  apex,  the  exposed  wing  is  small  relative  to  the  body,  so  that  the 
body  can  be  approximated  by  a  vertical  plane,  and  the  exposed  wing  ret  looted  in 
this  plane  forms  a  delta  wing.  This  'delta  wing'  may  now  be  thought  id  as  a 
fixed  wing  at  incidence.  Since  it  is  moving  downwards  with  velocity  -  r in 
dimensional  form  and  the  oncoming  flow  has  velocity  l'  ,  the  situation  is  equi¬ 
valent  to  a  wing,  in  a  flow  at  an  angle  of  incidence  i  =  -  r.  /F  ,  which  is  equal 
to  the  roll  rate  parameter  P  defined  in  the  previous  section.  Hence  an  avail¬ 
able  conical  flow  solution  will  give  the  unknowns  X(x^)  *”  the  cross-flow  plane 

x  =  x^  ,  say,  where  x^  is  sufficiently  small. 
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With  the  initial  solution  X(x  )  known,  we  use  a  marching  procedure  based 
on  the  Crank-Nicolson  method  for  advancing  the  solution  through  successive  cross- 
flow  planes  x^+ (  =  xm  +  hx  »  m  =  0,  i,  2,  ....  where  h^  is  the  step  length 
in  the  downstream  direction.  At  each  step  our  objective  is  to  satisfy  in  some 
cross-flow  plane,  the  stream  surface  and  pressure  conditions  (49)  and  (50)  at  the 
intermediate  points  on  the  sheet,  the  two  components  of  the  force  condition  (52), 
and  the  Kutta  condition  (53)  at  the  leading  edge,  giving  2n  +  3  equations  to 

solve  for  the  2n  +  3  unknown  components  of  X(x  ,)  .  These  equations  are 

-  m+ 1 

solved  by  an  iterative  method  described  in  detail  later.  The  method  starts  from 

a  guess  at  the  solution  vector  X(x  .)  and  uses  the  residuals  (ie  the  differ- 

-  m+ 1 


ences  between  the  left-  and  right-hand  sides)  of  equations  (49),  (50),  (52)  and 
(53)  to  improve  on  this  guess.  The  solution  at  x  =  x^  is  simply  the  conical 
solution  and  further  downstream  linear  extrapolation  from  the  two  previous  cross- 
flow  planes  is  used  to  generate  an  initial  guess. 

We  first  describe  how  the  quantities  appearing  in  (49),  (50),  (52)  and  (53) 

are  calculated  at  x  =  x  +  h  /2  from  the  known  solution  vector  X(x  )  and  a 

m  x  -  m 

guess  at  the  solution  vector  X(x  ,)  .  The  streamwise  derivatives  are  approxi- 

-  m+ 1 

mated  by  central-difference  formulae: 


-  f(xm  ♦  hx/2,  n) 
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For  the  remaining  terms,  linear  interpolation  is  used  to  provide  an  approximation 


to  X ( x  +  h  / 2)  : 
-  m  x 


X(x  +  h  12)  =  1[  X (x  )  +  X(x  )  1  +  OftT)  , 

-  m  x  m  -  m+1  1  x' 


(57) 


and  from  this  we  proceed  to  evaluate  the  quantities  needed  in  the  cross-flow 

plane  x  =  x  +  h  /2  . 

m  x 

As  indicated  earlier  we  use  an  intrinsic  coordinate  system  (q,  \Ji)  where  q 
is  defined  in  equation  (32),  to  represent  the  sheet  shape  in  eacli  cross-flow 
plane.  In  discretizing  the  sheet  itself  in  the  cross-flow  plane  we  require,  in 
particular,  an  accurate  representation  of  the  sheet  near  the  leading  edge  of  the 
wing  where  its  shape  changes  most  rapidly.  To  achieve  this  we  define  a  new 
variable  t  such  that 
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where  ?  is  the  quantity  introduced  in  section  3  and  X  is  a  constant  chosen  to 
give  the  required  angular  extent  of  the  sheet  at  the  initial  station.  In  the 
results  presented  in  the  next  section  X  is  chosen  so  that  at  the  initial  station 
the  angle  made  by  the  radius  vector  from  the  isolated  vortex  to  the  free  end  of 
the  sheet  with  the  real  axis  is  2  radians,  and  ■  is  given  by  (68). 

The  formula  for  the  sheet  coordinates  (35)  now  becomes 


Z  =  s  -  I  eM'f  dt  .  (59) 
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We  now  define  n  +  I  pivotal  points  to  be  equally  spaced  in  i  ,  with  step 

length  h  ,  say,  and  the  leading  edge  chosen  as  the  first  pivotal  point. 

Throughout  the  calculations  we  have  set  h  =0.1  .  Similarly  n  intermediate 

points  are  defined  with  the  same  spacing  in  t  but  commencing  at  the  point 

t  =  h  /2  .  The  representation  (58)  of  arc  length  was  introduced  by  Barsby""  ;  it 
T  2 

gives  n  ~  7,\t“/6  near  the  leading  edge  and  the  quadratic  dependence  upon  r 
ensures  an  adequate  density  for  the  pivotal  and  intermediate  points  close  to  the 
leading  edge.  Also  the  intermediate  points  in  successive  cross-flow  planes  are 
defined  by  constant  values  of  t  and  (58)  then  ensures  that  n  is  constant  for 
these  points  as  is  required  in  the  definition  of  the  streamwise  derivatives  on 
the  sheet. 


It  is  convenient  at  this  stage  to  adopt  a  numbering  system  for  the  quanti¬ 
ties  evaluated  at  the  pivotal  and  intermediate  points,  such  that  to  the  pivotal 
points  are  assigned  odd  numbered  subscripts  and  to  the  intermediate  points  even 


numbered  subscripts.  Thus,  for  example, 

intermediate  points  and  n„ .  ,  refers  to  n 

2  j- 1 


refers  to  q  evaluated  at  the  jth 
evaluated  at  the  jth  pivotal  point. 


From  X(x  +  h  /2)  we  know  the  values  of  y  at  the  j  t  h  intermediate  point, 
-  m  x 

i(>2j  ,  j  =  1  ( 1 )  n  ,  and  using  these  values  we  may  calculate  y.,._|,  j  =  l(l)n+l  , 
the  values  of  at  the  pivotal  points.  Because  the  sheet  must  leave  the  wing 


tangentially  we  have  t(’ j  =  0  and  the  remainder  are  found  using  four  point 
Lagrangian  interpolation  formulae.  Thus 


c 
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0 


3 


2j-» 


and 


2n-l 


2n+ 1 


(-  4 *  +  1 5i>.}  +  1  Oi)j^  -  *6)/20  , 

(-  * 2-  4  +  ^  +  9*  -  *,j+,)/l6  ,  j  =  3(1) n~ 1  , 

V  (60) 

(*2n-6  ~  5*2n-4  +  1  5*2n-2  +  5^2n)/16 

(“  5*2n-6  +  2 1  ^2n-4  "  35^2n-2  +  35*2„)/16  *  J 


The  sheet  coordinates  at  the  pivotal  points,  Z.  ,  ,  j  =  >(>)n+l  ,  may 

now  be  calculated  from  (59)  using  Simpson's  rule  to  achieve  comparable  accuracy 
with  the  interpolation  formulae  (60).  Thus 


and 


Z2j  +  1 


Z2.H 


h  e 

T 


2  j  —  I 


n- 


2j-l 


Mr,- 

+  4c  *”r) 


2j 


l^2 j  + 1 
+  o  J  n 


2j  +  l 


, j  =  1 ( 1 ) n  , (6 1 ) 


where  the  prime  denotes  a  derivative  with  respect  to  t  .  To  calculate  the 
coordinates  of  the  intermediate  points,  Z,;.,  j  =  KOn  ,  we  must  first  evaluate 


+ 


h  v’ 

i 

48 


*>•>  , 

Me  nj  +  19e  n, 


5e 


(62) 


This  formula  has  t  lie  same  order  of  accuracy  as  Simpson's  rule.  We  may  then  eval¬ 
uate  the  positions  of  the  remaining  intermediate  points  using  Simpson's  rule  in  a 
manner  similar  to  that  for  the  pivotal  points.  Hence 


h  ■ 

7=2  +  — — 

2j  2j-2  b 


r  i^2  j-1  ,  H'2j 

e  ‘h'j.2  +  4e  -J  n’ +  e 


11 2  j 


>  J 


!(l)n 


The  sheet  strength  is  specified  in  terms  of  the  parameters  t ■  ,  where 

«j  '  -  (If)  '  -  (iff),  j  ' 


The  values  at  the  pivotal  and  intermediate  points  are  related  by  equations 
like  (60)  for  the  angles  if  .  To  calculate  the  sheet  circulation  at  the  inter¬ 
mediate  points  we  have 
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..  I  I)A4>  , 

A<1  =  I  +  I  -  - - do  . 

J 

a 


To  evaluate  this  we  adopt  a  procedure  similar  to  that  which  gives  Z.,.  at  the 
intermediate  points.  Thus  we  first  evaluate  the  circulation  at  the  last  inter¬ 
mediate  point  using  an  integration  formula  analogous  to  (62),  and  for  the  other 
intermediate  points  we  add  on  the  contribution  from  the  rest  ot  the  integral 
which  is  calculated  using  Simpson's  rule.  The  circulation  at  the  intermediate 
points  is  then  given  by 


A4>„ 


2n 


48 


9«> 


+  | 

n+l  jn 


>K2n- 


g2n-2 


and 


A<i> 


h,  r 

2j  ~  *  T  f.’i 


i  i + 1  +  *\V2 


.  i  -  n- 1 i )  i 


The  sheet  coordinates  in  the  transtomud  plane  ■ .  ,  j  =  l(l)2n+l  ,  are 

calculated  from  the  corresponding  coordinates  in  the  c ross-l  low  plane  using 
equation  (4).  Hence,  in  non-d  (mens  i  ona  1  t  orm 


r, .  -  /.*  +  ►/*  -  i 

i  i  l 

where  7*  =  (Z. .  +  l/Z.)s/(l  *■  s')  .  (6  1) 

J  J  J 

Because  the  sheet  with  which  we  are  concerned  lies  entirely  in  the  1  irsl  quadrant, 
the  branch  of  the  square  root  with  positive  imaginary  part  is  used. 

To  calculate  the  complex  velocity  '.Hv'/a,-  at  any  point  in  the  trans¬ 

formed  plane  we  must  evaluate  the  integral  in  equation  (44)  numerically.  This 
must  be  considered  for  three  separate  cases: 

(i)  when  the  point  f,  is  not  on  the  sheet,  as  in  the  force  condition, 

Simpson's  rule  can  be  used  immediately  since  the  integrand  is  non-singular 
everywhere . 

(ii)  When  r.  =  1  ,  as  in  the  Kulta  condition,  the  integrand  is  singular  at 

u  =  0,  where  r,(n)  =  r,(o)  =  I  .  We  show  that  the  singularity  is  intolerable  and 
obtain  a  numerical  procedure  for  evaluating  the  integral.  We  note  lirst  that  the 
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last  two  terms  in  the  bracket  under  the  integral  sign  in  (54)  cancel  at  o  =  0 
when  4=1.  To  examine  the  behaviour  of  the  first  two  terms,  we  recall  that 
the  sheet  leaves  the  wing  tangential ly,  so  that,  for  points  on  a  sheet  close  to 
the  leading  edge 


Hence,  by  (4)  or  (63), 


s  +  o  , 


/ 2(s2  -  QoV 

\  s(s"  +  I  )  / 


4(o)  ~  C(o) 


Hence,  for  4  =  I  ,  the  first  two  terms  in  the  bracket  in  (54)  are 


-  4(o) 


4(o) 

4(o)  -  1 


/2s(s2  f  1)V 

\  (s2  -  Do/ 


Since  OA'l'/Do  is  finite  at  the  leading  edge,  by  the  Kutta  condition,  this  shows 
that  the  singularity  at  o  =  0  is  integrable.  To  remove  the  singularity,  we 
change  the  variable  of  integration  to  t  .  By  (58) 


o  -  7,U’t“76 


and  so 


I _>  d„  -  .  (M) 

\  c(o)  4(a)  -  1  /  y  3(S2  -  i)  ) 

With  t  as  variable  of  integration,  the  integrand  is  therefore  regular  and 
Simpson's  rule  can  again  be  used.  The  value  of  the  integrand  at  r  =  0  is  given 
by  combining  (64)  with  a  four-point  Lagrangian  extrapolation  formula  for  DA 4>/3o  : 

(FL  ■  £(’<!“).- ls^) ■  “» 


(iii)  When  DW/'Dr,  is  evaluated  at  f,  =  4^  on  the  vortex  sheet,  in  order  to 

calculate  q  and  q  from  (51),  the  integrand  is  singular  and  the  integral 
tm  n 

must  be  regarded  as  a  Cauchy  principal  value.  The  singular  behaviour  arises  from 
the  first  term  in  the  bracket  of  (54)  and  is  dealt  with  as  follows: 
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where  f  =  r. (o  )  and  the  partial  derivative  signs  merely  indicate  that  x  is 
constant.  The  first  term  on  the  right  is  regular  and  the  integral  in  the  second 
can  he  evaluated  analytical lv  as  a  Cauchy  principal  value: 


1 


Consideration  ol  the  simple  case  in  which  the  sheet  lies  along  the  real  axis  and 
I  -  '  •  .  tor  which  the  integral  should  he  real,  shows  that  the  imaginary 

part  of  the  logarithm  must  he  non-negative  and  less  than  2'i. 

Finally  in  the  numerical  formulation  of  the  problem  we  need  to  evaluate  the 
limit  in  equation  (52)  expressing  the  force  condition.  The  only  term  in  (54)  for 
■'W/V'.  which  gives  rise  to  a  singularity  in  3W/3Z  is  !'/2ni(,\  -  :  )  ,  so  we 
need  only  consider  t ho  limit  id 
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V  V  V 


as  Z  l  ends  1 Z^  .  expanding  -  r,  in  a  Taylor  series  gives 
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since 


i£ 

3Z 


1 


2s 


2 

+  s 


by  (4)  . 


This  completes  the  account  of  the  numerical  formulation  of  the  equations  gover¬ 
ning  the  model  and  we  now  proceed  to  describe  how  they  are  solved. 

We  have  explained  how  satisfying  the  discretized  boundary  conditions  at 

x  =  x  +  h  /2  gives  rise  to  2n  +  3  equations  in  the  values  taken  by  the 
m  x 

2n  +  3  unknowns,  d) .  ,  e.,  T  and  Z  ,  at  x  =  x  ,  .  We  have  denoted  the 
’  j  j  v  m+1 

vector  of  these  unknowns  by  X(x  ,)  .  When  a  trial  set  of  these  unknowns  is 

-  m+ 1 

introduced  into  the  set  of  2n  +  3  equations,  they  generate  2n  +  3  residuals, 
each  residual  being  the  difference  between  the  left-  and  right-hand  sides  of  an 
equation.  Let  us  denote  the  vector  formed  by  these  residuals  by  f(x(x  ))  • 

Our  object  is  to  modify  X  so  as  to  reduce  F  to  zero.  Note  that  neither  the 
ordering  of  the  unknowns  in  X  nor  that  of  the  equations  in  F  is  significant. 

r  .  r 

If  we  have  currently  a  vector  X  and  the  solution  is  X  ,  and  X  is 
close  to  X  in  some  sense,  we  can  use  the  multivariate  form  of  the  Taylor  expan¬ 
sion  to  write 

0  =  F (X)  =  F(Xr)  +  J  .  (X  -  Xr) 

where  J  is  the  Jacobian  matrix  of  the  system,  ie  the  matrix  with  SF^/SX^  in 
its  ith  row  and  j th  column.  If  the  sign  of  approximate  equality  is  replaced  by 
equality  and  the  resulting  equation  is  solved  for  X  ,  it  is  reasonable  to  expect 

that  the  solution  will  be  closer  (in  an  appropriate  sense)  to  the  true  solution 

r  .  .  .  r  + 1 

than  X  is,  so  it  is  taken  as  the  next  approximation  X  : 

Xr+1  =  Xr  -  J_1F(Xr)  .  (67) 


This  is  a  straightforward  multivariate  extension  of  the  Newton-Raphson  procedure 

for  the  iterative  solution  of  a  single  equation  in  one  unknown.  As  in  the  single- 

.  .  .  r  .  . 

variable  case,  convergence  is  quadratic  if  X  is  sufficiently  close  to  the 
solution  and  J  is  calculated  exactly.  The  success  of  this  method  however  does 
not  depend  upon  an  accurate  value  for  the  Jacobian  and  it  is  calculated  numeri¬ 
cally  from  the  initial  approximation  using  the  simple  forward  difference  formula 


Fi(X|,...,Xj_[,Xi  +  6,XM,...,X|<)  -  F  .  (X  , 


,X.,...,X.  ) 
_ j _ 


A. 


A 


♦ 


5  b 


where  k  =  2n  +  3  and  5  is  a  suitably  chosen  quantity  for  the  evaluation  of 

.  .  .  .  r+ I  r 

the  derivatives.  The  numerical  technique  used  for  calculating  (X  -  X  )  in 

(b7)  is  based  on  Gaussian  elimination.  Basically  one  wishes  to  solve  the  set  of 

linear  equations 

JX*  =  F  , 

,  .  r  r+ 1 

where  in  this  problem  X*  =  X  -  X  .  To  achieve  this  the  matrix  J  may  be 
written  in  the  form 

J  =  LI' 

where  L  and  U  are  Lower  and  upper  triangular  matrices  respectively  and  we 
have 

LUX*  =  F  . 

This  represents  the  two  triangular  systems 


LY  =  F  and  UX*  =  Y 


The  first  is  solved  for  Y  using  forward  substitution  and  t he  second  solved 

using  back  substitution.  This  method  for  solving  a  set  of  linear  equations  is 

29 

discussed  in  more  detail  m  Forsythe  and  Mo  1 e r  .  Once  the  solution  has  been 
obtained  for  X  ( xm+ 1 )  *  that  is  when  /?  ^  F  1 0  ^  ,  the  numerical  procedure 
continues  marching  downstream. 


When  using  the  solution  procedure  described  above  instabilities  can  be 
encountered  in  practice  which  manifest  themselves  as  oscillations  associated  with 
successive  downstream  solutions.  Ultimately  the  oscillations  lead  to  a  situation 
in  which  the  procedure  fails  to  converge.  Numerical  experiments  indicate  that 
the  oscillations  occur  when  the  errors  in  the  initial  approximation  which  is 
adopted  are  significant  in  relation  to  the  downstream  step  length.  In  the  present 
work  it  is  clear  that  the  conical  solution  adopted  as  the  starting  solution  is 
only  an  approximate  solution,  but  that  the  approximation  improves  as  x^ 
decreases.  However  as  decreases  it  is  found  that  in  order  to  maintain 

accuracy,  and  for  the  Newton  iteration  to  converge,  h  must  also  decrease  as 
the  scale  ef  the  vortex  system  is  diminishing.  As  a  consequence  it  is  difficult 
to  avoid  the  instability.  By  contrast  for  the  non-conical  flows  considered  in 
(17)  the  starting  solution  was  determined  exactly  as  a  conical  solution  and  the 
instability  did  not  occur. 


To  overcome  the  instability  in  the 
technique  is  employed.  This  consists  of 
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the  technique  described  above  and  then  calculating  X(Xp  +  h  /2)  from  (57).  This 

is  then  used  as  the  initial  solution  to  advance  the  solution  to  x_  +  3h  /2  and 

0  x 

an  averaged  solution  at  +  h^  is  found.  This  process  is  then  continued  down¬ 
stream.  This  smoothing  technique  damps  out  any  instability  very  effectively. 
However  it  has  only  the  accuracy  of  a  backward-difference  scheme  and  therefore  is 

less  accurate  by  a  factor  of  the  order  of  h  than  the  scheme  which  we  have 

x 

proposed  above.  We  therefore  choose  h^  =  0.01  ,  at  most,  when  the  smoothing 
technique  is  used  to  overcome  potential  instability,  whereas  much  larger  values 
of  h^_  may  be  used  when  the  original  scheme  is  employed.  In  this  way  the 
reduced  accuracy  of  the  smoothing  technique  is  accommodated.  The  smoothing  tech¬ 
nique  is  used  for  the  first  few  steps  until  it  proves  possible  to  obtain  a  con¬ 
verged  solution  for  a  larger  step  length,  or  until  the  errors  associated  with  the 
approximate  solution  have  been  damped  out.  Thereafter  the  original  full  central- 
difference  method  is  used. 

As  the  solution  is  advanced  downstream  the  number  of  iterations  required  at 
each  step  to  obtain  a  solution  decreases.  The  numerical  procedure  can  then  be 
speeded  up  considerably  by  only  recalculating  the  Jacobian  J  and  repeating  the 
triangular  decomposition  at  the  start  of  the  Newton  iteration  when  the  number  of 
iterations  required  at  the  previous  step  is  greater  than  some  specified  number. 
Typically  it  was  recalculated  when  the  number  of  iterations  needed  at  each  step 
was  greater  than  6  and  it  was  by  no  means  unusual,  especially  for  large  distances 
downstream,  for  the  procedure  to  march  for  many  steps  before  recalculating  J  . 

6  RESULTS 

6. 1  Parameter  values 

We  wish  to  compare  results  obtained  by  the  present  method,  using  slender- 

body  theory  for  a  wing-body  combination,  with  Harvey's  measurements^'  on  a  delta 

wing.  (The  present  method  is,  of  course,  equally  applicable  to  wings  of  any 

slender  planform.)  Regarding  the  delta  wing  within  the  context  of  slender-body 

theory,  in  which  there  is  no  upstream  influence,  we  note  that  the  only  length- 

scale  is  the  local  semi-span  s^(x)  •  If*  therefore,  for  two  delta  wings  of  the 

same  apex  angle,  each  rolling  at  zero  incidence,  there  are  cross-flow  planes  for 

which  the  values  of  S2s,(x)/U  are  the  same,  the  flow  fields  in  these  planes  will 

d 

be  the  same,  when  scaled  in  terms  of  the  length  s^(x)  and  free-stream  speed 

U  .  For  the  wing-body  combination,  this  similarity  no  longer  applies,  since  the 
body  radius  r  provides  a  further  length-scale. 

To  make  a  comparison  with  the  results  for  the  wing  alone,  the  present 
method  should  be  applied  for  a  vanishingly  small  value  of  r  .  In  terms  of  the 
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non-dimensional  variables  introduced  in  section  4,  this  means  an  indefinitely 
large  value  of  s  (the  semi-span  divided  by  r  )  and  a  vanishingly  small  value  of 
P  (=  -  ilr/U),  such  that  the  product  Ps  remains  finite  and  equal  to  -  iJs  /U  .  In 
practice,  we  are  limited  to  finite  values  of  s  and  P  ,  so  calculations  have 
been  made  for  Ps  equal  to  the  values  of  -  Sls^/U  appropriate  to  the  measure¬ 
ments,  for  several  small  values  of  P  .  It  is  apparent  from  these  that  the  solu¬ 
tions  depend  little  on  P  ,  unless  Ps  is  small,  and  some  features  of  the 
solutions  can  be  extrapolated  to  P  =  0  with  confidence. 

The  values  of  P  chosen  are  0.04,  0.02  and  0.01.  The  semi-apex  angle  of 
Harvey's  wing  is  10°.  As  explained  in  the  previous  section,  the  downstream 
integration  starts  with  an  initial  solution  appropriate  to  a  delta  wing  at  an 
angle  of  incidence,  a  ,  equal  to  P  .  The  values  of  the  incidence  parameter 
a/tan  y  corresponding  to  these  three  rates  of  roll  are  therefore  0.227,  0.113 

and  0.057.  The  last  two  of  these  are  below  the  range  for  which  solutions  could 

.  4  5 

be  found  by  Smith  ,  using  a  polar  representation  of  the  sheet,  but  Barsby 

obtained  solutions  using  the  same  intrinsic  coordinate  system  as  the  present 
treatment.  At  the  lowest  value  of  a/tan  y  the  sheet  has  a  marked  inflexion, 
and  at  only  slightly  lower  values  Barsby  found  that  he  could  only  obtain  solu¬ 
tions  by  allowing  the  separation  point  to  move  inboard  from  the  leading  edge.  We 

note  in  passing  that  the  present  results  cannot  be  generalized  to  other  values  of 

22 

y  ,  in  the  way  that  the  results  of  Hamn  and  Mishne  ,  using  an  isolated  vortex 
model  for  the  rolling  delta  wing,  can  be  generalized. 

The  number,  n  ,  of  points  representing  the  sheet  was  chosen  to  be  12.  A 

,  .  .  4 

similar  degree  of  detail  was  found  adequate  by  Smith  and  the  need  to  carry  the 

calculations  to  large  values  of  s  in  reasonable  time  argued  against  an 

unnecessarily  large  value. 

The  marching  procedure  in  all  these  cases  was  started  in  thi  cross-flow 
plane  for  which  s  -  1.1  ,  that  is  x  =  x^  =  0.1/y  and  the  quantity  •  U'od  in 
the  definition  of  q  ,  (58),  was  specified  as 

7  =  (s  -  I)  +  0.01(x  -  •  1 b8) 

In  the  immediate  neighbourhood  of  the  initial  station,  l  grows  like  (x  -  1)  as 
required  for  a  conical  approximation,  and  increases  quadrat ical ly  as  we  march 
further  downstream. 

In  his  experimental  investigations  Harvey  obtained  detailed  pressure  meas¬ 
urements  on  the  wing  surface  at  three  locations  on  the  wing,  namely  x/c  =  0.333, 
0.55(3  and  0.778,  where  c  is  the  root  chord  of  the  wing,  and  at  the  three 
different  roll  rates 


0.1,  0.2  and  0.3  . 
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'•sd(c) 

U 


Thus  we  have  results  for  nine  values  of  Ps  available  from  the  experiments 
namely  0.0333,  0.0556,  0.0666,  0.0778,  0.0999,  0.1112,  0.1556,  0.1668  and  0.2334. 
These  are  the  values  for  which  a  detailed  analysis  of  the  results  obtained  has 
been  carried  out. 

6 . 2  Vortex  sheet  shapes  and  circulations 

In  Fig  4  we  show  the  vortex  sheet  shape  and  vortex  position,  scaled  with 
the  local  semi-span  s  ,  for  P  =  0.01  at  various  stieamwi se  stations.  These 
streamwise  stations  correspond  with  the  values  taken  by  Ps  in  Harvey's  experi¬ 
ments.  The  vortex  system  grows  downstream  in  a  non-conical  fashion,  as  can  be 
seen,  but  the  form  specified  for  f(x)  in  (68)  ensures  that  the  angular  extent 
of  the  sheet  remains  almost  constant.  This  appears  to  be  a  happy  coincidence, 
for  although  the  quadratic  behaviour  of  £(x)  may  reflect  some  underlying  struc¬ 
ture  to  the  solutions  we  have  been  unable  to  expose  it.  This  constant  angular 
extent  is  also  maintained  in  the  results  for  P  =  0.02  and  0.04  ,  and  in  Fig  5 
the  sheet  shapes  for  the  three  rates  of  roll  are  shown  for  certain  of  the  values 
of  Ps  .  These  demonstrate  that  the  different  body  size  has  little  influence 
upon  the  sheet  shapes  once  the  semi-span  is  large  compared  to  the  body  diameter, 
and  in  particular,  that  P  =  0.01  is  small  enough  for  the  effects  of  the  body  on 
the  sheet  shape  to  be  small  for  these  values  of  Ps  .  The  positions  of  the  iso¬ 
lated  vortices  are  again  shown  in  Fig  6  and  are  seen  to  remain  close  together  for 
the  different  values  of  P  .  Thus  the  vortex  position  does  not  depend  signific¬ 
antly  on  r/s  .  Harvey  measured  the  total  head  distribution  in  the  planes 
x/c  =  0.5,  0.6  and  0.7  for  the  same  angular  velocities  as  in  his  pressure  measure¬ 
ments.  The  smallest  of  the  contours  of  total  head  surrounding  the  vortex  cores 
which  he  was  able  to  draw  are  also  shown  in  Fig  6,  for  the  highest  angular 
velocity.  These  show  that  vortex  cores  in  the  experiment  were  slightly  lower 
than  the  isolated  vortices  in  the  calculation  and  about  5a  of  the  semi-span 
further  inboard.  This  discrepancy  is  about  the  same  as  that  found  in  the  case  of 
a  delta  wing  at  incidence,  weich  is  provisionally  ascribed  to  the  effects  of 
secondary  separation  in  Ref  4. 

2  2. 

Hanin  and  Mishne"  have  carried  out  a  ealeulalion  for  the  separated  tlow 
past  a  rolling  delta  wing  using  an  isolated  vortex  model,  following  Brown  and 
Michael".  The  resulting  vortex  positions  are  also  shown  in  fig  6.  As  tor  the 
ease  of  a  fixed  delta  wing  at  incidence  the  predicted  in  .lit  of  the  vortex  agrees 
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well  with  that  given  by  the  vortex  sheet  model  while  the  lateral  distance  of  the 
vortex  from  the  wing  centre  line  is  significantly  overestimated. 

The  total  circulation  (non-dimensional)  on  the  vortex  sheet  and  isolated 

2  .  ... 
vortex,  divided  by  Ps  ,  is  plotted  in  Fig  7.  In  dimensional  variables  this 

2 

corresponds  to  the  total  circulation  divided  by  -  12s  .  The  graph  shows  this 

quantity  rising  fairly  quickly,  as  Ps  increases,  to  a  plateau  of  roughly  con- 

2 

stant  height,  showing  that  the  circulation  is  about  1.5  Ps  for  s  >  5  .  Hanin 
2'5 

and  Mishne~“  found  that  the  first  term  in  the  expansion  (for  small  rate  of  roll) 

2 

for  the  circulation  of  the  isolated  vortex  is  proportional  to  -  fts  ,  but  it  is 
not  clear  how  significant  this  correspondence  is. 

6 . 3  Wing  surface  pressures 

From  Bernoulli's  equation  (14)  and  relation  (19)  the  pressure  coefficient 
is  given  by 

C„  =  -  2<J>  -  f<l>2  +  4>2]  +  2pif>  z  -  $  y  >  . 

P  x  y  z>  (  y  zJ  ) 


The  velocity  components  and  <J>  may  be  calculated  directly  from  the  real 

and  imaginary  parts  of  3W/3Z  but  must  be  evaluated  using  a  difference 

formula  for  the  downstream  derivative.  In  order,  therefore,  to  determine  the 

pressure  coefficient  at  a  point  (x^  +  h^/2,  y,  0)  on  the  wing  we  evaluate  the 

vector  X(x  +  h  /2)  from  X(x  )  and  X(x  +  h  )  by  linear  interpolation  as  in 
-  m  x  -  m  -  m  x 

(57)  and  hence  determine  <S  and  <t>  .  To  determine  the  downstream  derivative 

y  z 

of  the  potential,  $  ,  in  the  cross-flow  plane  x  +  h  /2  ,  we  calculate  the 

x  m  x 

potential  v  =  R{W)  for  the  same  value  of  y  in  the  planes  x  =  x 

m 

x  =  x  +  h  using  the  result  obtained  by  integrating  (54),  namely 
m  x 


The  constant  of  integration  in  (69),  which  could  in  principle  depend  on  x  ,  has 
been  chosen  to  be  zero,  so  as  to  ensure  that  ,  and  therefore  C  ,  tends  to 

X  P 

zero  at  inf inity . 


c 
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Wo  then  calculate  T  using  the  central  difference  formula  (56).  When 

x 

evaluating  (69)  care  must  be  taken  to  ensure  that  the  correct  branch  of  the  com¬ 
plex  logarithm  is  taken.  This  is  determined  by  the  criterion  that  the  logarithm 
is  continuous  everywhere  except  across  the  vortex  sheet  where  it  jumps  di scout in- 
uously  by  2  ?i  i  .  The  numerical  techniques  used  to  evaluate  the  integral  in  (69) 
are  similar  to  those  described  in  section  5. 

Detailed  pressure  distributions  over  the  wing  surface  for  1’  =  0.04,  0.02 
and  0.01  are  shown  in  Figs  8  to  13,  in  which  the  pressure  coefficient  is  plotted 
against  the  distance  from  the  centre  line,  non-dimensional ised  with  respect  to 
the  overall  semi-span.  The  pressure  coefficient  was  not  calculated  over  the  body 
as  it  makes  no  contribution  to  the  rolling  moment,  and  so  an  indication  of  the 
relative  size  of  the  body  and  wing  is  given  by  the  ordinate  at  which  each  pressure 
distribution  terminates.  Each  figure  i r  for  a  particular  value  of  Ps  and  so  the 
smaller  body  sizes  correspond  to  the  smaller  values  of  P  .  For  the  lower  values 
of  Ps  the  suction  peaks  on  the  upper  surface  are  very  high  and  narrow;  as  Ps 
increases  the  height  of  these  peaks  falls  slightly,  but  they  broaden  considerably. 
These  peaks  also  illustrate  that  the  solutions  for  different  values  of  P  but  the 

same  value  of  Ps  are  very  close.  For  the  higher  values  of  Ps  the  experimen- 

2  1 

tally  measured  pressure  coefficients  of  Harvey  are  superposed.  In  most  of  these 

cases  the  agreement  is  good  only  away  from  the  suction  peak  and,  as  for  the  case 

of  a  wing  at  incidence,  in  the  region  of  the  peak  the  experimental  curves  are 

lower  and  broader.  For  the  last  case,  however,  Ps  =  0.2334  ,  the  agreement  is 

2  ’ 

very  good.  The  pressure  distributions  calculated  by  Hanin  and  Mishin*"’''  are  not 
shown  as  their  calculated  suction  peaks  are  considerably  higher  and  narrower  than 
those  shown  here. 

Several  factors  may  be  responsible  for  the  differences  between  the  calcula¬ 
ted  and  measured  pressure  distributions.  It  could  be  that  the  presence  of  the 
body  is  affecting  the  calculated  values.  This  effect  would  he  largest  for  the 
smallest  value  of  Ps  ,  vet  Fig  9  shows  that  the  effects  ot  a  tour-fold  change  in 
body  size  are  negligible  compared  with  the  difference  between  the  calculated  and 
measured  pressures  at  Ps  -  0.08  .  At  the  streamwise  station  nearest  to  the 
trailing  edge,  x/c  =  0.7/8  ,  the  upstream  influence  of  the  trailing  edge,  neglec¬ 
ted  in  the  theory,  could  be  expected  to  affect  the  measured  values.  Results  for 
this  station  are  shown  in  Figs  9,  II  and  13,  and  it  is  difficult  to  identity  any 
consistent  discrepancy  in  these  three  figures.  It  appears  theretore  that  upstream 
influence  of  the  trailing  edge  is  fairly  small  on  this  wing  at  this:  station,  which 
is  about  1.25  semi-spans  ahead  of  the  trailing  edge. 
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The  flow  on  the  upper  surface  underneath  the  vortex  is  directed  outboard 
relative  to  the  vortex  and  encounters  a  rising  pressure  between  the  vortex  and 
the  leading  edge.  This  pressure  rise  will  tend  to  produce  a  secondary  separation 
of  the  boundary  layer,  not  represented  in  the  theoretical  model.  Certainly  the 
very  steep  gradients  and  large  pressure  rises  predicted  by  the  theory  for  the 
smaller  values  of  Ps  seem  incompatible  with  attached  boundary  layers,  either  lam¬ 
inar  or  turbulent.  On  the  other  hand,  Harvey  suggests  that  no  secondary  separation 
occurred  in  his  experiments,  on  the  grounds  that  appreciable  pressure  rises  were 
observed  between  the  vortex  and  the  leading  edge.  He  attributes  the  absence  of 
separation  to  a  thinning  of  the  boundary  layer,  produced  by  the  rotation  of  the 
wing.  At  the  time  of  Harvey's  work,  almost  all  the  available  pressure  distribu¬ 
tions  on  wings  with  leading-edge  vortices  had  been  measured  in  circumstances  in 
which  the  boundary  layer  was  laminar  at  secondary  separation,  and  in  these 
circumstances  little  pressure  rise  is  observed.  Since  then,  appreciable 
recompression  of  the  flow  outboard  of  the  vortex  has  been  found  on  wings  at 
incidence,  in  cases  where  the  boundary  layer  is  turbulent  in  this  region,  see  for 
example  Ref  30.  Typically,  the  turbulent  boundary  layer  sustains  a  larger  pressure 

rise  before  separating,  and  then  separates  further  outboard,  than  the  laminar 

.  .  .  30 

layer.  Comparison  with  typical  shapes  of  measured  suction  peaks  suggests  that 

those  in  Figs  II  and  13  and,  probably  Fig  9  are  associated  with  turbulent  second¬ 
ary  separation  and  that  those  in  Fig  10  and,  possibly.  Fig  12  correspond  to 
laminar  secondary  separation.  This  would  imply  that  transition  occurred  between 
x/c  =  0.556  and  0.778,  it-:  at  a  Reynolds  number  of  between  5  and  7  *  10“*  based  on 

distance  from  the  apex.  This  is  somewhat  lower  than  the  corresponding  Reynolds 

30 

number  for  wings  of  aspect  ratio  1  at  moderate  incidence  .  In  the  absence  of 

direct  observations,  it  is  thought  likely  therefore  that  secondary  separation  was 

present  in  Harvey's  experiments.  In  the  more  thoroughly  examined  flow  over  a 

4 

wing  at  incidence,  it  has  been  argued  that  the  discrepancy  between  the  calculated 
and  measured  suction  peak  is  due  to  secondary  separation,  largely  on  the  basis 
that  transition  has  a  marked  effect  on  the  form  of  the  suction  peak.  We  therefore 
ascribe  most  of  the  difference  between  the  calculated  and  measured  suction  peaks 
to  the  absence  of  a  representation  of  secondary  separation  in  the  theoretical 
model . 

6.4  The  local  rolling  moment 

The  quantity  of  most  practical  interest  is  the  rolling  moment  induced  by  the 
rolling  motion,  the  damping  in  roll.  The  present  approach  lends  itself  most  rea¬ 
dily  to  the  calculation  of  the  local  rolling  moment,  ic  the  rolling  moment  acting 
on  a  narrow  element  of  the  wing  between  two  neighbouring  cross-flow  planes, 


oiy 
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divided  by  the  distance  between  them.  This  is  also  the  quantity  derived  by 
Harvey-*  from  his  pressure  measurements.  The  local  rolling  moment  can  be  expres¬ 


sed  conveniently  in  coefficient  form  by  referring  it  to  the  product  of  the  kinetic 

2  ...... 

pressure  (jpU-)  and  the  square  of  the  local  span.  This  coefficient  is  divided  by 

the  roll  rate  parameter,  -  i2s/U  ,  to  give  a  quantity  which  varies  slowly  with 

roll  rate,  and  would  be  the  constant  stability  derivative  in  a  linear  model  of  the 

flow.  Thus  we  have  a  derivative,  C  ,  of  the  local  rolling  moment  coefficient 

ip 

given  by 


rolling  moment/unit  chord 
1  U  - 


2  2 
ipU  4s 


where  the  quantities  are  dimensional.  Hence 


ccP  =  ik  j  ACP(X>  y)  s  d@ 

0 

where  the  quantities  are  now  non-dimens i onal . 

For  the  case  of  a  rolling  slender  delta  wing  in  attached  flow  this  may 

2  1 

easily  be  evaluated  (on  Harvey  )  to  give 


a 

C  =  -  r  tan 

>P 


The  at tached-f low  rolling  moment  for  the  wing-body  combination  was  evaluated 
numerically,  using  a  subroutine  QA04A  ,  from  tie  Harwell  Subroutine  Library, 
which  automatically  takes  care  of  the  square  root  singularity  at  the  leading  edge. 
These  results  agree  very  closely  with  the  formula 


0 


P 


tan  y  (s'4  -  1) 

ii  4 

s 


tan 


s  + 


s(s-  +  I) 


(70) 


which  is  derived  by  differentiating  the  overall  rolling  moment  coefficient  calcul- 

■>  )  ... 
ated  by  Lomax  and  Heaslet  .  This  coefficient  is  sketched  against  Ps  in  Fig  14 

for  attached  flow  and  for  P  =  0.01  .  For  the  rolling  wing-hody  combination  the 

roe f f i c i on t  qnie.ly  rises  1 rora  zero,  overshoots  the  value  for  the  delta  wing  and 

then  settles  down  very  close  to  this  line. 


Simple  flow  models,  the  isolated  vortex  model  (Hanin  and  Mishne“)  and  the 
leading-edge  suction  analogy  model  (see  the  Appendix)  may  also  be  employed  to 
give  values  for  the  local  rolling  moment  coefficient  for  the  delta  wing.  These 
are  also  sketched  in  Fig  14.  To  determine  the  effect  of  the  body  upon  the  solu¬ 
tions,  calculations  were  carried  out  using  the  isolated  vortex  model  applied  to  a 
rolling  wing-body  combination.  In  these  calculations  the  numerical  method  dis¬ 
cussed  earlier  in  section  5  was  employed,  except  that  no  vortex  sheets  were  pres¬ 
ent  and  therefore  only  the  Kutta  condition  and  zero  force  condition  were  used. 

These  results,  for  P  =  0.01  and  0.02  are  shown  in  Fig  14.  Comparing  the  results 

22 

with  those  of  Hanin  and  Mishne  for  P  =  0  we  see  that  at  the  larger  values  of 
Ps  both  sets  of  solutions,  with  and  without  the  body,  are  in  close  agreement 
although  they  differ  considerably  for  the  lower  values  of  Ps  .  To  attempt  to 
eliminate  the  effect  of  the  body  we  may  linearly  extrapolate  the  results  of  the 
isolated  vortex  model  for  P  =  0.02  and  0.01  to  P  =  0  ,  corresponding  to  r  =  0  , 
the  rolling  delta  wing.  In  this  way  we  do  indeed  reproduce  the  results  of  Hanin 
and  Mishne,  except  for  Ps  =  0.033  where  the  effect  of  the  body  is  significant. 

This  may  be  easily  seen  in  Fig  14,  since,  for  Ps  0.04  ,  the  difference  between 

the  curves  for  P  =  0.02  and  0.01  is  equal  to  the  difference  of  the  results  tor 
P  =  0.01  and  those  of  Hanin  and  Mishne  for  P  =  0  .  We  should  therefore  expect  a 
reasonable  approximation  for  a  rolling  delta  wing  to  be  given  by  extrapolating  the 
calculated  results  to  P  =  0  .  Similarly  an  analysis  may  also  be  carried  out 
using  the  leading-edge  suction  analogy  and  the  results  are  sketched  in  Fig  14. 

The  details  of  the  method  are  given  in  the  Appendix.  Those  again  show  that 
extrapolation  of  the  curves  for  P  =  0.02  and  0.01  to  P  =  0  reproduce  the 
results  for  the  case  when  no  body  is  present.  For  the  full  vortex  sheet  model, 
therefore,  we  should  expect  a  reasonable  approximation  for  a  rolling  delta  wing  to 
be  given  by  extrapolating  the  results  for  low  values  of  P  to  P  =  0  . 

This  local  rolling  moment  coefficient  was  evaluated  for  the  full  vortex 
sheet  model  for  a  range  of  values  of  Ps  ,  with  P  =  0.01,  0.02  and  0.04  .  The 

integration  technique  used  employed  Simpson's  rule  over  100  and  200  integration 

points,  together  with  Richardson  extrapolation  to  increase  the  accuracy.  As  a 
check  upon  the  accuracy  some  of  the  integrations  were  repeated  using  200  and  400 
integration  points  with  no  appreciable  difference  in  the  results.  The  rolling 
moment  coefficients  are  shown  in  Fig  15.  For  fixed  Ps  the  suction  peaks  are 
always  slightly  higher  for  the  larger  values  of  P  .  Also,  since  they  are  close 
to  the  Leading  edge  the  fact  that  they  are  higher  has  a  significant  effect  on  t  lie 
rolling  moment.  This  effect  far  outweighs  the  extra  contributions  from  the  larger 
wing  length  for  the  smaller  values  of  P  ,  corresponding  to  the  smaller  body  sizes. 
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Consequently  the  rolling  moment  coefficients  are  in  general  larger  for  large  P  , 
except  in  those  cases  for  which  the  wing  span  is  not  much  bigger  than  the  body 
diameter. 

As  we  mentioned  earlier,  we  should  expect  a  reasonable  approximation  for  a 
rolling  delta  wing  to  be  given  by  extrapolating  the  results  for  the  wing-body 
combination  to  P  =  0  .  This  extrapolated  curve  is  indicated  in  Fig  15. 

In  Fig  16  the  experimental  results  of  Harvey  are  given,  together  with  the 
extrapolated  results  for  P  =  0  from  the  different  theories.  We  see  that  the 
isolated  vortex  model  considerably  overpredicts  the  magnitude  of  the  local  rolling 
moment  coefficient.  This  is  because  the  calculated  suction  peaks  arc  too  far  out¬ 
board  and  too  high.  The  leading-edge  suction  analogy,  however,  gives  better 
agreement.  As  can  be  expected  from  the  pressure  distribution  shown  in  Figs  8  to 
13,  the  results  for  the  vortex  sheet  model  agree  well  with  the  experimental  value 
for  i’s  =  0.2334  but  not  so  well  for  smaller  values  of  Ps  .  Any  difference 
between  the  theoretical  and  experimental  pressure  distributions  close  to  the  lead¬ 
ing  edge  of  the  wing  will  of  course  be  magnified  when  calculating  the  rolling 
moment.  When  secondary  separation  is  present  the  modification  to  the  pressure 
distribution  in  flows  of  this  type  is  such  that  the  normal  force  is,  in  general, 
predicted  correctly.  When  this  is  the  case  the  rolling  moment  cannot  be  predicted 
correctly.  It  is  not  surprising  therefore  that  the  closest  agreement  between 
theory  and  experiment  occurs  for  t lie  experimental  results  associated  with  the 
station  closest  to  the  trailing  edge  and  the  highest  rate  of  roll.  For,  in  spite 
of  the  upstream  effect  of  the  trailing  edge  in  this  case,  the  fact  that  the 
separation  is  turbulent  in  the  experiments  leads  to  a  pressure  distribution 
similar  to  that  predicted  by  the  theoretical  calculations. 

7  CONCLUSIONS 

The  vortex  sheet  model  of  leading-edge  separation  has  been  successfully 
adapted  to  calculate  the  positions  and  strengths  of  the  leading-edge  vortices,  the 
pressure  distribution  over  the  wing  and  the  local  rolling  moment  of  slender  delta 
wings  mounted  on  cylindrical  bodies  in  steady  rolling  motion  at  zero  incidence. 

Comparisons  with  the  experimental  observations  of  Harvey  suggest  that  the 
model  describes  this  flow  as  well  as  it  describes  the  flow  past  a  slender  wing  at 
incidence.  However,  because  the  discrepancies  occur  close  to  the  leading  edge, 

1  he  prediction  of  rolling  moment  due  to  roll  is  less  accurate  than  the  prediction 
of  lift  due  to  incidence.  The  vortices  in  both  flows  are  found  to  lie  inboard  of 
their  calculated  positions  to  1  to  57  of  the  local  semi-span;  and  the  form  of  the 
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measured  pressure  distribution  agrees  best  with  the  calculation  when  the  Reynolds 
number  is  higher  and  the  boundary  layer  can  be  assumed  turbulent  at  secondary 
separation.  These  are  also  the  conditions  in  which  the  prediction  of  the  rolling 
moment  is  most  accurate. 

The  results  confirm  the  measurements  and  the  predictions  of  simpler 
22  31 

models  *  that  leading-edge  separation  increases  the  roll  damping.  The  present 
predictions  are  more  realistic  than  those  of  the  isolated  vortex  model.  They  also 
show  that  the  presence  of  the  body  (in  a  configuration  of  fixed  total  span) 
increases  the  roll  damping  unless  the  exposed  wing  span  is  very  small. 
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Appendix 

APPLICATION  OF  THE  LEADING-EDGE  SUCTION  ANALOGY 

The  results  presented  earlier  for  the  rolling  moment  coefficient  have  also 

been  compared  with  some  results  obtained  using  an  empirical  method  derived  by 
19 

Polhamus  for  calculating  the  lift  on  a  slender  wing  and  later  extended  to  the 

31. 

rolling  case  by  Boyden  .  In  this  theory  the  nonlinear  lift  force  acting  on  the 
wing  is  represented  by  the  force  arising  from  the  singularity  at  the  leading  edge 
of  the  wing.  In  his  paper  Boyden  used  lifting  surface  theory  to  obtain  the 
linear  and  nonlinear  lift  contributions  but  in  the  present  work  we  just  use  the 
results  predicted  by  slender  body  theory,  for  consistency  with  the  other  models 
used . 

32 

From  Robinson  and  Laurmann  the  leading-edge  suction  torce  per  unit 
chord  is  given  by 

> 

F  =  real(-  ,vC“)  (A-l) 


where  C  is  defined  by 


fA 

dz 


f(X  -  s)  +  - 

v7.  -  s 


in  the  neighbourhood  of  the  leading  edge,  Z  =  s  .  In  the  present  case,  from  (4) 
we  have,  in  dimensional  notation 


s2  -  r2  | 


DZ 


+  0(1) 


2s  (s"  +  r")  iX  -  s 


and  hence 


dWA  3W 

_ A  _  _ A 

ax  ~  a.-. 


=  1) 


2  2 
s "  -  r“ 


+  0(1) 


2s  (s“  +  r")  ►/.  -  s 


where  dW^/Or,  is  given  by  (11).  Hence 
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2  2 

s'"  -  r* 

■>  T- 
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(-*.  =  I) 


(A-2) 
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and  so,  through  (A-l)  the  suction  force  F  acting  on  the  leading  edge  may  be 
calculated.  This  force  acts  in  a  direction  tangential  to  the  wing  surface  and, 
according  to  Polhamus,  the  nonlinear  force  acting  on  the  wing  due  to  the  leading- 
edge  separation  is  given  by  rotating  the  force  F  through  90°  so  that  it  acts 
perpendicular  to  the  wing.  Assuming  that  this  force  still  acts  through  the  lead¬ 
ing  edge  the  magnitude  of  the  local  moment  per  unit  chord  is  given  by  s | F  j  . 

Now 


rolling  mnmpnf/nm'  r  nhnrH 


where  the  quantities  are  dimensional .  Hence  from  (A~l)  we  obtain  the  rolling 
moment  coefficient  due  to  the  suction  on  both  sides  of  the  wing 


(C  ) 

Ip  nonlinear 


|Re(C2) 

Ps2 


(A—  3) 


where  the  quantities  are  now  non-dimensional.  C  is  given  from  (A-2)  as 


C 


1) 


J _ 

+  I) 


(A-4) 


and  so,  through  (A-l)  the  force  F  acting  on  the  singularity  may  be  calculated. 

The  total  rolling  moment  coefficient  is  then  found  by  adding  on  the  contri¬ 
bution  from  the  attached  flow  potential.  Now,  at  the  leading  edge,  from  (13),  we 
have 


ac 


C=1  IT  COS  6 


— ^ —  ( tt  20  j  +  sin  20  () 


where  cos  0  ^  «  2s / ( 1  +  s  )  .  Since 


0|  =  cos 


•1  2s 


1  +  s 


~  -1  tt 

=  2  tan  s  -  —  , 


we  obtain  from  (A-3) ,  (A-4)  and  (A-5) 


(A-5) 
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LIST  OF  SYMBOLS 

unit  vectors  in  tangent  plane  to  the  vortex  sheet 

constant  (see  Appendix) 

local  rolling  moment  coefficient 

vector  of  residuals 

=  -  3A$/3t| 

n 

downstream  step  length 

step  length  in  cross-flow  plane 

unit  vectors 

function  of  x  used  to  maintain  a  constant  angular  extent  for  the 
sheet 

number  of  intermediate  points  on  the  sheet 

unit  outward  normal  to  the  wing-body  combination  and  vortex  sheet 
pressure 

roll  rate  parameter  (=  -  Qr/U) 

velocity  vector  relative  to  the  wing-body  combination 

relative  velocity  tangential  to  the  sheet  in  the  cross-flow  plane 

relative  velocity  normal  to  the  sheet  in  the  cross-flow  plane 

radius  of  body 

wing-body  semi-span 

time 

velocity  vector 

complex  potential 

attached  flow  complex  potential 

vector  of  unknowns  for  Newton  iteration 

Cartesian  coordinates 

=  y  +  iz 

position  of  isolated  vortex  in  cross-flow  plane 
wing  semi-apex-angle 

difference  operator  across  the  vortex  sheet 
coordinate  on  the  vortex  sheet 

coordinate  in  complex  transformed  cross-flow  plane 

position  of  isolated  vortex  in  the  transformed  plane 

=  cos  '  (2rs/(r^  +  s^)) 

scalar  quantities  (section  3) 

coordinate  system  (section  3) 

density 


P 
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LIST  OF  SYMBOLS  (concluded) 

0  arc  length 

vortex  sheet  (section  3) 
t  stretched  coordinate  on  sheet 

l1  velocity  potential 

<Jj  angle  between  tangent  to  sheet  and  wing 

The  same  symbols  have  been  used  for  dimensional  quantities  before  section  4 
for  quantities  non-dimensionalized  using  r  and  U  later  in  the  main  text, 
the  Appendix  both  dimensional  and  non-dimensional  quantities  are  used. 
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